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report describes the selection, sizing, and optimization of energy storage/power 
ining systems (ES/PC) for the electromagnetic propulsion systems. An ES/PC is the 
e between an electromagnetic thruster and a primary power source. The ES/PC ac- 
le relatively low, continuous power output from the primary power source, and con- 
to a series of high power pulses and delivers it to the thruster. The ES/PC 
s the primary power requirements while maximizing the thruster performance. 

lection, sizing, and optimization can be achieved if an optimization criterion 
ana if the ES/PC is examined within the context of the integrated propulsion sys- 
seiected mass minimization as the optimization criterion. To develop the selec- 
hoaology, we identified the requirements imposed on the ES/PC and parametrically 
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:d) ; The propulsion systems considered covered a range of missions, 

thrusters, primary power options, and ES/PC options. The missions 
ationkeeping and orbit transfer from lower earth orbit (LEO) to geo- 
(GEO) of large payloads (up to 25 metric tons). The thrusters considered 
inductive thruster (TIP), and the electric rail gun (ERG). The primary 
:: solar cell arrays, solar thermophotovoltaic, nuclear thermonic, nuclear 

iclear Brayton cycle, nuclear RanJcine cycle, and nuclear magnetohydro- 
le ES/PC options were: capacitive, normal inductive, superconducting 
ertial based systems. 


he selection method over the range of missions, thrusters 
found that for each mission/^nruster combination, there 
lizes the mass of the propulsion system. 
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Mp^p/(exp(AV/v^) - 1) 


(20) 


In order to solve this eqnatioi.. the mass of the propellant and the dry 
engine mass must be related to the thrnster output parameters as described 
by Equations 17 to 19. 

The second equation derived from the m.ision requirements is that the 
total thrust time, t^, is related to the mission time, t^ by: 


(21) 


where a = the thruster duty cycle. 

The thruster duty cycle is defined by the mission, and Equation 21 uniquely 
defines the total thrust time. 

No further independent relations among the thruster output parameters 
can be derived from the mission. The two remaining equations must be de¬ 
rived from propulsion system optimization. 

2 3.2 Propulsion System Optimization 

Optimized design of electromagnetic propulsion systems involves the 
minimization of deployed mass. The criterion of mass minimization applied 
to system selection and sizing is illustrated as follows. 

Mathematically, minimization of the total propulsion system mass may be 
thought of as minimizing a function (the propulsion system mass) of several 
independent variables (the propulsion system outputs) subject to certain 
anziliary restrictions (the mission requirements) The propulsion system 
mass is the sum of the subsystem masses'- 

«tot = Mpp - ^ (22) 

which can be expressed as a function of the propulsion system outputs since 
all the subsystem masses are uniquely related to these variables (Equations 
17 to 19). Then Mtot written as: 

“tot = (23) 
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If these relations were known, inspection of Equations S, 8 and 12 shows 
that all the major propulsion subsystem masses (and other performance re¬ 
quirements) are uniquely related to the propulsion system outputs as: 


“PP= 


(17) 




(18) 


(19) 

prp sys prp sys p d 

Therefore, the performance characteristics of the major subsystems 
must be known before an integrated propulsion system can be selected, sized, 
or designed in an optimal manner. In order to illustrate this more clearly 
we next examine the generic problem of propulsion system selection and 
sizing. 


2.3 PROPULSION SYSTEM SELECTION AND SIZING 

Tbe selection and sizing of an optimal prop nl si on system for a specific 
mission is a complex task. The basic principles are simple. Propulsion 
system selection and sizing may be thought of as the process of choosing a 
combination of m^, f and t^j which meet the mission objectives and 

satisfy the optimization criteria. Also, the process of choosing the pro¬ 
pulsion system outputs results in a unique definition of the size and per¬ 
formance of each of the propulsion subsystems. 

Mathematically, the process may be thought of as simultaneously solving 
four equations with four unknowns (i.e.. the propulsion system outputs). In 
order to find unique solutions for these output parameters, four independent 
equations which relate the parameters to one another must be developed. The 
process by which this is accomplished is illustrated in the following 
paragraphs. 


2.3.1 The Mission 

The mission provides us with two of the four required equations. The 
first equation is derived from the momentum change which must be imparted to 
the payload mass. Momentum must also be imparted to the dry engine mass and 
the propellant. The complete equation is: 
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(11) 


“tank = kjMp^p 

and the propellant system mass can be written as: 

“prp sys ~ ^®p^^d 

where k 2 ^ seen from equation 12 , thruster characteristics do 

not directly affect propellant mass requirements. 

2.2.4 The Thruster 

The thruster converts the electrical energy delivered to it from the 
ES/PC into kinetic energy of the pulse mass. The thruster is required to do 
this at a frequency and for a total time consistent with the mission 
requirements. The manner in which the thruster accomplishes this task has a 
major impact on the other subsystems in the propulsion system. Thruster 
performance has a first order effect on the ES/PC and primary power systems. 
The efficiency of the thruster will determine the energy storage and average 
power requirements of the system. The plysics of thruster operation will 
determine the current and voltage levels and pulse duration required of the 
ES/PC. 

Dnder 'steady' average thrusting conditions the thruster produces iden¬ 
tical pulses and its performance should be largely independent of the pulse 
frequency or the total operating time. This implies that thruster per¬ 
formance is only dependent on pulse mass and exhaust velocity. The modeling 
relationships which are required may then be expressed as 


’'th* nth^-p'^e^ 
Jth = Jth<“p'^e> 
'^th = 

tp = tpdBp.V^) 


(13) 

(14) 


11 






“espc = M,,pc(Ep.J.V.f.t 

where J = ES/PC output current 
V = ES/PC output voltage 
f > operating frequency 
tp s pulse duration 
t^ * thrusting time 

The efficiency of the ES/PC system will also depend on these factors and may 
be expressed as: 

%spc = nespc(Ep'J'V.f.tp,t^) (9) 

The relationships described by Equations 8 and 9 have not been extensively 
investigated. The determination of these relationships was one of the 
primary goals of this prograia. 

Examination of Equations 8 and 9 irdicates that four of the six param¬ 
eters on which ES/PC performance depends are determined by thruster perfor¬ 
mance. These four parameters are Ep, j, V and tp. Thruster modeling must, 
therefore, not only address pulse energy requirements (or efficiency) but 
must also treat current, voltage, and pulse duration. 


The primary functions of the propellant system are to store the propel¬ 
lant and to provide the reaction mass for the propulsion system. The 
propellant mass is a significant mass component of the total propulsion 
system and in some missions will dominate propulsion system mass. The 
propellant mass is related to propulsion system outputs by: 


Mp,p - 

The propellant storage container mass may be a considerable fraction 
of the propellant system mass, depending on the propellant type (gas or 
solid) and on the way it is stored (i.e., gas propellant may be stored 
cryogeni cal ly or in high pressure tanks). In general, the container mass is 
expressed as a function of the propellant mass as: 
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Ibe range of mission parameters nsed in this stndy is simimarized in 
Table 3. 


TABLE 3 

SUMMARY OF MISSION CHARACTERISTICS 



Payload Mass 

AV 

Mi ssi on 

Duty 


(kg) 

(m/s) 

Time 

Cycle 

Stationkeeping 

2000-25000 

50-185/yr. 

1-10 yrs 

0.1-1 

Orbit Transfer 

2000-25000 

5800 

30-300 days 

1 


2.2 INTEGRATED ELECTROMAGNETIC PROPCLSION SYSTEM 

A generic electromagnetic propulsion system is illustrated in Figure 2. 
It consists of four major subsystems or mass components: a primary power 

source, an ES/PC, a propellant system (propellant and its storage con¬ 
tainers) and the thruster itself. All of these subsystems are required to 
produce the thruster output, a pulse of propellant mass, m^, at an average 
exhaust velocity, v^. The propulsion system produces these mass pulses 
repetitively at frequency, f, for a total time, t^. 

The parameters mp, v^, f and t^j are the propulsion system outputs. 
Therefore, system performance must be modeled with these parameters as the 
independent variables. 



(594) 


Figure 2. Generic Electromagnetic Propulsion System. 

















A/M for a generic class of space structures which can be deployed in a 
single shuttle launch is presented in Table 1 and ranges from 0.08 to 0.598 
m^/kg. The AV requirements calculated for this generic class of space 
structures are presented in Figure 1. 


TABLE 1 

CHARACTERISTICS OF A GENERIC O^ASS OF SPACE STRDCTDRES^ 


SINGLE SHUTTLE 


LADNCH 




MASS 

AREA/MASS 

EFFECTIVE 

GENERIC GLASS 

SIZE 


(kg) 

(m2/kg) 

A/M(m^/kg) 

I PLATE 

MEDIUM 

(100 

m) 

1618 

4.014 

0.201 

STHDCTDRE 
W/0 BLANKET 

LARGE 

(250 

m) 

3672 

11.055 

0.553 

II PLATE 

SMALL 

( 30 

m) 

1334 

0.43.> 

0.438 

STHDCTDRE 

MEDIUM 

(100 

m) 

11350 

0.572 

0.572 

W/BLANKET 

LARGE 

(150 

m) 

24420 

0.598 

0.598 

III MODDLAR 

SMALL 

( 15 

m) 

2300 

0.165 

0.091 

ANTENNA 

ICDIDM 

( 60 

m) 

8375 

0.433 

0.113 


LARGE 

(200 

m) 

18017 

1.980 

0.236 

IV SERIES OF 

SMALL 

(2) 


7500 

0.826 

0.085 

ANTENNAS 

MEDIUM 

(3) 


11250 

0.802 

0.085 


LARGE 

(4) 


15000 

0.764 

0.084 


The total AV requirements for all stationkeeping functions range from 
50 to 184 m/s/yr mission times ranging from 1 to 10 years. 


2.1.2 Orbit Transfer 

The AV for orbit transfer is 5800 m/s. The AV required for LEO to GEO 
transfer was investigated in a previous program and the results are pre¬ 
sented in Table 2. The mission times of interest range from the minimum 


^Smith, W. W. and Clark, J. P., 'Study of Electrical and Chemical Propulsion 
Systems for Auxiliary Propulsion of Large Space Systems,' NASA CR-165502, 
1981. 









The AV required to correct the payload inclination caused by solar and 


lunar perturbations is given by: 


AV/year = A0 AVq£q (2) 

where A0 = required payload orbital plane inclination change (rad/yr) 

A0 A sin e cos e 

A » constant whose value depends on the body creating the effect 
(0.74 for solar and 1.615 for lunar) 

6 = inclination of equatorial plane with respect to the ecliptic 

or lunar orbital planes 

^^GEO “ tangential velocity of GEO payload (3074.7 m/s) 

Substituting back into Equation 2 we get: 


AV/year » -(3074.7/(360/2n))[(0.74 sin 0 cos 9)^^^^^ 

+ (1.615 sin e cos e)ian*r^ (3) 

The AV ranges froai 41.2 to 51.8 n/s/yr. 

The AV required to correct eccentricity changes caused by solar radia¬ 
tion effects is given by: 


AV/year = (3 skit/2A) (ajt/2 sin (oji/2))0/sin“^ p) (4) 

wh<*re s =■ solar constant at 1 AO » 4.5xl0‘'^kg/(iii-s^) 
k - (l-fo)(A/M) 

a « average reflectivity of payload - 0.3 

A/M » area-to-nass ratio (m^/kg) 

K a mean angular velocity of earth motion around sun 
= 1.99xl0“^rad/s 
a » duty cycle 

8 ^ eccentricity ratio 

* ■§-, e* a maximum allowable eccentricity 

Cp a peak eccentricity that would occur if initial orbit 
was circular and no stationkeeping corrections were 
applied. 

8 defines the stationkeeping accuracy required and is set to zero to calcu¬ 
late the maximum AV. AV is a function of a and A/M% a may be as high as 1. 





SECTION 2 


ES/PC ANALYSES APPROACH 


The ES/PC will be examined within the context of an integrated propnl 
sion system. The process of propulsion system selection begins with identi 
fication of the mission. The propulsion system is then sized to accompl is: 
that mission in the most economical manner (the propnlsion system i 
optimized for the mission). Optimized design can, in principle, always bi 
accomplished if criteria for optimization are established. We assnmed tha' 
the optimization criterion is minimization of the integrated propnlsioi 
system mass. 


2.1 MISSION CHARACTERISTICS 

All space propnl sion systems are required to accomplish the same gen¬ 
eric task: imparting a velocity increment, AV, to a payload mass, xi 

some mission time, tg,* These mission parameters are the basic, independeni 
variables which constrain propulsion system selection and sizing, and t 
which propulsion system performance must be related. The range of missioi 
parameters of interest must be identified. 


2.1.1 Stationkeeping 

On-orbit stationkeeping of geosynchronous payloads is required to re¬ 
move accumulated disturbances caused by: longitude drift caused by triaxi- 
ality of the earth, long term increase in inclination of the payload orbital 
plane due to luuar and solar perturbations, and change in eccentricity 
caused by solar radiation effects. Stationkeeping is accomplished by sup¬ 
plying external torques which are provided by the propulsion systeai. Th< 
AV imposed on the propulsion system by these disturbances were calculated. 

The AV due to longitudinal drift is expressed as: 

AV/year = (1.72 sin (2k+15)) (1) 

where k is the longitude in degrees East of Greenwich. The AV range is froe 
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SECTION 1 


INTEODOCTION 

An energy storage/power conditioning system (ES/PC) is the interface 
between an electromagnetic thrnster and a primary power source in a space 
propulsion system. An ES/PC accepts tbe relatively low, continuous power 
output from a primary source, converts it to a series of high power pulses 
and delivers it to the thrusters. Electromagnetic thrusters require peak 
power levels well above the average power level required to accomplish 
mission objectives. The ES/PC is required to minimize the primary power 
source requirements while maximizing thruster performance. 

Development of electromagnetic propulsion systems has continued for 
over twenty years. Host of the work has concentrated on the thrusters. As 
thruster technology matures, more emphasis must be placed on the ES/PC and 
on system integration. This report describes the results of a program 
conducted at lAP Research, Inc. to establish guidelines by which program 
planners can optimally select and evaluate the performance of integrated 
propulsion systems. The guidelines allow the program planner to make 
tradeoffs among the propulsion system components and compare alternative 
concepts. 

The propulsion systems considered in this program covered a range of 
missions, thrusters, and primary power options. The missions were: ten 

year stationkeeping and orbit transfer from low earth orbit (LEO) to geosyn¬ 
chronous orbit (GEO) of large payloads (up to 25 metric ton). The thrusters 
were: the Teflon pulsed plasma (TPP), the magne topi asmadynamic (MPD), the 

pulsed inductive (PIT), and the electric rail gun (ERG). The primary power 
options were: solar cell arrays, solar thermophotovol taic, nuclear ther¬ 
mionic, nuclear thermoelectric, nucl ear Rankine cycle, nucl ear Bray ton cycle 
and nuclear magnetohydrodynamic (MHO). 

Selection, sizing, and optimization of an ES/PC can be achieved if an 
optimization criterion is established and if the ES/PC is examined within 
the context of an integrated propulsion system. To develop the selection 
guidelines, we identified the requirements imposed on the ES/PC and paramet¬ 
rically related the ES/PC performance to the propulsion system outputs. 






SECTION 3 


INTBGEATED FKOPDLSION SYSTEM mFORMANCE AND MODELING 

In this section we describe the characteristics of the subsystems in 
the integrated propulsion system and present the performance and mass models 
of these subsystems. 

3.1 PRIMARY POWER CHAEACIERISTICS 

The primary power source was assumed to contain all the components 
necessary to deliver electrical power, at some characteristic voltage and 
current level, to the ES/PC The performance characteristics of several 
power sources were studied and mass models were developed. Tie characteris¬ 
tics and mass models are described in detail in Appendix A. A snamary is 
presented in this section. 

The characteristics of interest are the output constraints (current and 
voltage), the minimum and maximum power levels available, and the mass. The 
characteristics were determined by conducting a literature survey and con¬ 
tacting technical experts to confirm the results and make refinements where 
needed. The literature contained considerable point design data but very 
little parametric or limit data. The primary power mast is best described 
as: 

“pp ■ -.p^pp >“> 

where m^p i* the specific mass. The specific mass was modeled by curve¬ 
fitting the data gathered in the literature search for power levels ranging 
from 2 hW to 7 MV. The results are presented in Table 4. All power sources 
investigated in this study can be theoretically designed for any power 
level. If a power level lower than the minimum specified in Table 4 is 
required, the mass is set to the minimum. For high power levels 0500 hV) 
multiple shuttle launches may be required. These effects are shown in 
Figure 4. 

The output voltage range of these systems is presented in Figure 5. 
There is a tradeoff between ohmic losses in the busses and insulation degra¬ 
dation due to the high voltage. The optimum output voltage cannot be deter¬ 
mined until the mission is defined. 
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TABLE 4 


SPECIFIC MASS OF THE PRIMARY POWER SYSTEMS 


POWER SOURCE 

SPEaFIC HASS 

<“sp) 

MODEL* MINIHIM MASS 




(kg) 

Solar Cell Arrays 

0.013 + 40.61/P„„. 

PP 

>1000 w 

PP 

54 

S ol ar 

TlieraiopliotoTol taic 

0.005 + 4.6xl0*‘^P 

PP 

1000 

W < Ppp < 10.000 W 

6 

Nncl ear 

Thermoelectric 

0.021 + 370.3/P„„. 

PP 

Ppp 

> 1000 W 

3 92 

Nad ear 

Thermionic 

0.018 + 514/Ppp. Ppp > 10,000 W 

700 

Nuclear Brayton 
Cycle 

0.019 + 190.6/P_. 

PP 

^‘pp > 

1500 W 

220 

Nuclear Ranhine 
Cycle 

0.013 * 178.6/P^ 
PP 

Ppp 

> 1000 w 

191 

Nu^^ear MHD 

10-4 + 1990.1/P„ 
PP 

^’pp 

> 100,000 W 

2000 
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Figure 5. Power Output Voltage Range of Primary Power Systems. 


3.2 IHRDSTER CHARACTERISTICS 

Hie performance characteristics of the four electronagnetic thrusters 
were identified and their inputs were modeled. Hie characteristics and 
models are described in detail in Appendix B. A simimary of the characteris¬ 
tics and a brief description of the operation of each thruster is presented 
in the following paragraphs. 


3.2.1 Pulsed Inductive Thruster (PIT) 

The PIT consists of a flat, spiral coil through which an energy supply 
is discharged. A cross-sectional view is shown in Figure 6. A thin layer 
of propellant gas is injected over the face of the coil. When this propel¬ 
lant layer has achieved the desired density, the energy supply is discharged 
through the coil. 


9 , 
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Figure 6. Cross-sectional View of Pulsed Inductive Tiraster (PIT), 


in the opposite direction to current flow in the coil, and the two repel 
each other. The sheet accelerates away froa the coil in the axial direction 
and produces the desired iapul se. The acceleration of the propellant ceases 
at some distance fron the primary coil, when the magnitude of the primary 
magnetic field is insufficient to induce current flow in the propellant. 


3.2.2 Teflon Pulsed Plasma Thruster (TPP) 

The TPP thruster is illustrated in Figure 7, and consists of a set of 
parallel electrodes separated hy several Teflon fuel bars. The Teflon 
propellant is held against the electrodes by a constant force spring. A 
solid state ignitor plug is located in the cathode of the thruster, between 
the Teflon bars, and an energy supply is connected directly to the thruster 
electrodes. 

In operation, the energy supply is charged and supplies s potential to 
the thruster electrodes. This applied voltage is sustained by the thruster 
until a pulse is supplied to the ignitor plug initiating a discharge between 
the electrodes. 'Cie mini-discharge from the ignitor plug produces enough 
electrical conductivity between the thruster electrodes for the main energy 
supply to discharge. An arc forms between the electrodes, ablates and 
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Figure 7. Teflon Poised Plasma Thruster (TPP). 

ionizes material from the face of the Teflon fuel bars. The ionized propel¬ 
lant is then accelerated by gas dynamic and electromagnetic forces produced 
from the interaction of the arc current and its self-generated magnetic 
field. After the energy supply has discharged, the cycle can be repeated. 

3.2.3 Magnetoplasmadynamic Thruster (MPD) 

The HPO thruster has a coaxial cylindrical geometry as shown in 
'igure 8. It consists of an outer cylindrical anode and an axial cathode. 
The back wall of the thruster chamber is an insulating backplate which 
contains propellant injection holes. A power supply is connected across the 
electrodes. Propellant is injected into the chamber through the backplate 
.t a constant flow rate. The presence of the propellant gas between the 
ilectrodes causes a breakdown, forming a plasma. The high radial plasma 
;urrent (lO^A) induces a toroidal self-generated magnetic field which accel- 











3.2.5 Tirnster Performance Modelins 
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Figure 9. Electric Rail Gun (ERG). 


sconplisbed if a large body of data existed (either experinental or analyt- 
:al) for maxiBtai efficiency thrusters and if this data contained values 
ix the input requireaents over a wide range of and v^> A siaple nul ti- 
■riate curve fitting technique could then be used to develop the required 
slationships for E^, J. V and tp> This was the approach used for aodel ing 
le electric rail gun.^ This siaple date-based approach conld not be ap- 
Lied to the other thrusters, because not enough data was available over a 
de range of a^^ and v^. An alternate approach was developed based on 

ting perfoxBan( 



















































Tlie efficiency of the thruster is defined ss the kinetic energy divided 
by the pulse energy. For the PIT end TPP. the efficiency is: 

n - 1/(1 + 2 e./y2 + 2k/(mJ/3^^)) (32) 

Ihe efficiency for the MPO is: 

T, - 1/(1 f le^/v^ + 2ki/{iJ/2v^) + (33) 

Ihe efficiency for the ER(j^ is: 

n - 13.67 (34) 

3.2.5.2 Current/Pul se Tine 

An sversge current csn be defined ss: 

= Ul^ J^dt/tp]^/2 (35) 

which csn be further reduced by using the electroBSgne ti c inpulse equstion: 

apVp - (L'/j^dt)/2, (36) 

where L' is the srisl inductsnce grsdient for the device. Conbining Equs- 
tions 35 snd 36 we obtsin: 

“ (2 «pV^/(L'tp))l/2 (37) 

Equstion 37 nust be reduced further to obtsin the nodel of interest (s func¬ 
tion of and only). Therefore L' snd t^ hsve to be expressed ss 

functions of n^ snd v^. This is discussed in Appendix E The resulting 
sversge current snd pulse durstions sre 

tp - 13 np^''^/We (38) 

snd 

J.vg = 148.25 (39) 


for the PIT. 
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For the TPP, the expressions xre: 




(40) 


snd 


J^yg = 1.07 x10^«P^/3 (41) 

For the MPD, the pnl se dors tios (tp) is sn independent parsneter since the 
MPD is operated in a qnasi-steady mode. The MPD cnrrent is: 

= 2.42x10^ 

For the ERG^> the expressions are: 

tp = 4.45xl0“^mpl/3^^ (43) 

and 

J^yg - 2.2xJ0‘*mpl/3 (44) 


3.2 .5 .3 Voltage 

The thrusters considered in this study^ are all electromagnetic and are 
'current driven' devices (i.e., the acceleration forces are related to cnr¬ 
rent, not voltage). These devices do not require a high voltage to accel¬ 
erate the pulse mass (although some require a relatively high voltage for 
the formation of the current carrying plasma and to produce the required 
short pulses). The voltage is, therefore, a derived parameter in electro¬ 
magnetic thrusters and is indicative of the impedance of the device. 

The average voltage, V^yg> i* computed from the pulse energy, the pulse 
time, and the average cnrrent and can be written as: 

'^svg - 

Substituting Ep, t^ and J^yg *»y corresponding equation for each 

thruster we obtain the following. 
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For the PIT, 


- 2.59x10"‘*»p^^3v^^ + 4.61il04mpl^3 + 7.52xl0-2v^ (46) 

For the TPP, 


V^g . 0.311 ■*■ 7.5 x10^«P^^2/v^+87.2 (47) 


For the MPD, 

'^evg “ 2.07xl0"^ip^^^Ve^^^ + 1.16xl0^npl/^v^~^^2+o 495 v^l/2.^28.3 (48) 

For the EHG, 


V 


*vg 


0.032 » 

P « 


(49) 


3.2.5.4 Snaaxry 

Thruster chsxscteristics sxid perforaance were related to the propulsion 
system outputs. The models deweloped for the thrusters (except for the EKG) 
were based on an understanding of the physical principles of operation and a 
number of assumptions concerning the performance of optimized thrusters 
combined with awailable performance data. The models are extremely versa¬ 
tile, permitting rapid evaluation of thruster input requirements over a wide 
range of output conditions. 

3.3 PROPELLANT SYSTEM (BARACTERISTICS 

The propellant system mass is expressed as 

“prp sys “ K Vd '12) 

where E is a constant greater than one and accounts for the storage contain¬ 
ers and propellant feed system. The storage container mass varies over a 
vide range, depending on the propellant type (gas/solid) and the storage 
media (gas propellant may be stored cry oge ni ca 1 1 y or in high pressure 
tanks). Thus, the exact storage container mass cannot be determined until 
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The input constraints and output requirements of the ES/PC are essen¬ 
tial in identifying the concepts and components capable of performing the 
required functions. The constraints and requirements are presented in 
Table 6. These were computed from the primary power output and thruster 
input requirements. More than one mass, size, and efficiency model may be 
required for the ES/PC concepts and components, because the parameter values 
span such a wide range (up to 5 order of magnitude). 

TABLE 6 

ES/PC INPUT CONSTRAINTS AND OUTPUT REQOIREilENTS 



PARAMETER 

VALUE 

Output 

Ep (hJ) 

0.1 - 400 


V (kV) 

0.06 - 27.5 


I (kA) 

15.0 - 1170 


tp Us) 

1.0 - 5800 


f (Hz) 

0.04 - 19 

Input 

Ppp (kW) 

0.2 - 7800 


V (V) 

28 - 1000 


The energy store in an ES/PC is the main component of the system. The 
input conditioners are selected to charge it efficiently and the output 
conditioners are selected to deliver the stored energy in the required form 
i.e., to shape the pulse). Three types of energy stores were considered; 
capacitive, inductive, and inertial. 

Three capacitive-based ES/PC concepts were considered. The first con¬ 
cept consisted of a simple capacitor with a closing switch as the output 







assumed tiie total compulsator mass to be 2.5 times tiie rotor mass, which can 
be expressed as; 


“comp = 10 Ep/(av2) (69, 

The compulsator size can be derived from the rotor size. It will have 
the same length and twice the radius. The rotor radius, r, is related to 
the rotor tip speed by: 


= 2«fmr (70) 

Combining Equations 65, 67, and 70 we obtain: 

r = 8Pvtp/„ (71) 

The length (b) of the rotor can be calculated from the rotor mass and radius 
and is given by: 




(72) 


These equations can be used to calculate the mass and size of a compulsator 
optimized for a set of independent parameters (tp, Ep, and f) specified by 

the mission. During selection and sizing of a compulsator, it should be 

Itept in mind that; 

(1) a short pulse requirement may lead to an impr acti cal ly small rotor 
radius. A practical minima pulse duration is 100 ps to 500 qs. 

(2) a long pulse, on the other hand, may lead to a short rotor (thin 
disk) which is dynamically unstable. In this case, the tip speed 
must be decreased until the 1 ength-to-di ame ter ratio is satis¬ 
factory (say 1) and the mass and radius of the rotor must be 
recalculated. 

(3) Finally, surface current density must not exceed the maximum 

allowable limit of 10 MA/m in order to control the armature wind¬ 

ing bond shear stress. If the current density exceeds this value, 
the length of the rotor must be increased and the mass and rotor 
radius must be recalculated. 

We exercised this methodology over the range of independent parameters 
of interest and some of the results are presented in Figures 22 and 23. The 


42 







The pulse width is related to the electrical frequency, of the 

alternator by: 

tp = l/(l6f^) (65) 

because the compulsator coaipresses the half cycle pulse width of the alter¬ 
nator by a factor of 8 (aastning el ectronagnetic thrusters act resistively). 
The electrical frequency is related to rotational speed by: 

fe = P/2 (66) 

where f^ i* the mechanical rotor speed and P is the number of field poles. 
The factor of 2 can be removed from Equation 66 if either polarity can be 
used. The electrical frequency becomes: 

f, = f„ P (67) 

The minimum pulse width can be decreased by increasing either the rotor 
speed or the ntmiber of poles or both. The maximum rotor speed is stress 
limited to 150 m/s. The maximum number of poles is limited by the air gap 
distance between the rotor and stator. As the pole-to-pole separation 
approaches the air gap distance, the excitation field leakage increases. At 
some point, the leakage becomes so great that the addition of poles is 
futile. 

The mass components of a compulsator are: the rotor (the energy 
store), the back iron (return path of the magnetic flux) and the field coils 
(generate the magnetic field). The mass of the rotor is given by: 

“rot ' 

where v is the rotor tip speed. The energy stored is calculated from the 
ratio of pulsed to stored energy (a) and is limited to less than 5 percent 
by peak power constraints. The back iron has the same cross-sectional area 
as the rotor and its mass should be approximately equal to the rotor (both 
are ferromagnetic and have similar densities). The field coils are sized to 
deliver the required output voltage and their mass is relatively small. We 






(2048) 


Figure 21. Compalsator Cross Sectioa.^^ 

We condncted s literstore survey s&d developed s methodology to cslcu- 
iste the size snd mass of s compulsator. The methodology and the limita¬ 
tions of the compulsator are discussed below. 

The fundamental limitations to the compulsator were investigated by 
Weldon et al.^'^ The limits can be divided into three groups: those dealing 
with the effects of load characteristics, those limiting output power, and 
those limiting minimimi pulse width. 

The load affects the ability of the compulsator to compress flux and 
produce short pulses. The compulsator was found to reduce the basic alterna¬ 
tor half cycle pulse width by a factor of 8 for a resistive load and by a 
factor of 4 for a capacitive load. 

The peah output power is limited by several factors. The thickness and 
dielectric strength of the air gap between the rotor and stator windings 
have a significant effect on peak power. The air gap must be minimized in 
order to minimize the inductance and thereby maximize output current. This 
reduces the voltage breakdown strength of the gap (and the maximna permis¬ 
sible output voltage of the machine). A second limitation on the output 
power is the shear strength of the insulation system used to bond the stator 
windings to the stator and rotor windings to the rotor. The interaction 
between the compulsator discharge current and the excitation field causes a 
tangential force on the conductors. This force slows the rotor, converting 
inertial energy to electrical energy, and results in a tangential shear 
stress on the insulation bond. The maximm strength of the insulation was 
found to be 28 MPa (4 ksi) which corresponds to a maximum surface current 
density of 10 MA/m. 


40 






“ref ' 60 (64) 

where = refrigeration power. 

This aoes not include the mass of radiators required to reject the heat. 

3 .4 .3 .4 Ihe Compulsator 

The compensated pulsed alternator (compulsator) is an inertial energy 
store. It combines both the energy store and the output conditioning func¬ 
tions in one device. It provides high-voltage, high-current pulses by 
utilizing the principles of magnetic induction and flux compression. The 
compulsator is a single-phase alternator with a stationary field and a 
rotating armature. The armature winding and an identical stationary (or 
compensating) winding are connected in seriea. At at least one point per 
cycle the inductance of the armature circuit is minimized. Ihe variable 
armature inductance leads to a flux compression action. When coupled with 
alternator action this produces high current pulses into the load. The 
basic operation is illustrated in Figure 20. A simplified cross-sectional 
view of a compulsator is presented in Figure 21. 



Figure 20. Compensated Pulsed Alternator ( Comp nl sat .^r) 


^^Weldon, W. F.. Bird. W. L.. Driga, M. D., Tol k, K. M.. Rylander, H. G.. 
and Woodson, E E, 'Fundamental Limitations and Design Considerations for 
Compensated Pulsed Alternators', (£M University of Texas at Austin. 
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auili&ry eqaipnent. Ilie analysis and characterization of superconducting 
systems is described in detail in Appendix D. 

Soper conducting systems have four major components, the inductor, the 
current leads, the dewar. and the refrigerator. Ihe mass of the inductor is 
determined by the choice of geometry, the materials, the magnetic field, and 
the structural support system. When suitable choices are made for these 
parameters the specific energy of the inductor is found to be: 

Wind “ 3.33 X 103j/kg (61) 

The current leads carry current from the normal enviroment into the 
cryogenic environment of the inductor. Ihe current leads are a major mass 

component. The size of current leads is determined by a tradeoff between 

lead mass and resistive losses. The lead mass is then related only to the 
current carried and is given by: 

“lead * ^00 I (62) 

The dewar provides the cryogenic environment for the inductor and is 
designed to limit heat flow from ambient to cryogenic. The mass of the 
dewar is detemined by the size of the inductor and by the loads it must 
carry. The mass of the dewar and insulation is given by: 

“dew ‘ Adew («3) 

where A^ew * surface area of the dewar. 

The remaining component in a superconducting system is the refrigera¬ 
tor. ‘Die long duration of propul sion missions requires that a closed cycle 
refrigeration system be used. The size of the refrigeration system is 
primarily dependent on the heat losses into the cryogenic system. The 
losses result from the 'downleads.' the dewar. and the pulse losses in the 
superconductor. This energy must be removed by the refrigeration system at 
cryogenic temperature and rejected at ambient temperature. The mass of the 
refrigerator is primarily related to the refrigeration power required at the 
cryogen temperature. The mass is: 
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Hie characteristic diaeasions of a capacitor caa be calculated from the 
Bass. Hie mass of the dielectric is half the capacitor mass. Hie deasity 
aad Yolme of the dielectric are kaowa. Assmiag the capacitor is a cube, 
the characteristic dimeasioa mil 1 be the cube root of the volnme. 

Hie dielectric thickaess caa be calculated from Eqnatioas S3. 56, 57. 
The miaimiai possible dielectric thickaess was assumed to be 2.5 aicroas.^^ 
If the calculated dielectric thickaess is less thaa the miaiaum specified, 
thea it is set to the miaimum aad the mass of the capacitor recalculated. 

Hie efficieacy of capacitors was assiaed to be 90%. that is 10% of the 
eaergy stared is dissipated ia the teminatioas aad bus bars. 

3.4.3.2 Normal laductors 

laductors are used ia pulse power systems either as eaergy stores or as 
pul se-f oxmiag cimipoaeats. Whea used as eaergy stores, laductors are attrac' 
tire because they prow ide a coastaat curreat source characteristic. Hiis 
characteristic is yery iaportaat for some electromagnetic thrusters. The 
disadyantage of inductors is that they are 'lossy' aad inefficient. When 
used as pul se-f oiBing coaponents, inductors control the shape aad duration 
of pulses deliyered from capacitors. Hiey can be used either ia a normal 
'LC combination or in a multi-staged pulse-forming network. As pulse- 
forming components, inductors carry current for a short time aad therefore 
the resistiye losses inherent ia inductors, are not so important. 

Sizing inductors for either function (energy storage or pulse forming) 
is complex and must account for theoial, mechanical, and geometrical con¬ 
straints. Hie effects of these constraints on the size, mass, and effici¬ 
ency of normal inductors are discussed in detail in Appendix C. 

We are primarily interested in minimum-mass inductors. Hie minimum 
mass configuration is found by simultaneous solution of the equations de¬ 
scribing the constraints (thermal, mechanical, and geoetrical). A closed 
form solution is not possible and a numerical approach must be used. A 
computer program was deyeloped to find the solution. The inductance, cur¬ 
rent, and conduction time were treated as the independent yariables aad 
input to the computer program. The material properties such as resistiyity, 
density, and maximum stress must also be supplied. The outputs of the 
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(E/E^) 


(57) 


Substituting Equation 57 into Equation 55, we get the specific energy 
of dielectric for voltage reversal greater tliaa 20 percent: 

w - e^(lt/p)(N.e0*0^®*^/2.5)~®-^^^E§/2 

and for voltage reversal less tban 20 percent: (58) 

w - 8o(k/p)N^-^^^E^/2 

To maxinize tbe specific energy (ainiaize tixe aass). Equation 58 shows 
that IcE^/p should be aaxiaua. An evaluation of the ccaaonly used dielectric 
aaterials in capacitors is presented in Table 8. The aost attractive 
dielectric is polypropylene film. 


TABLE 8 

DIELECiaiC PROPERTIES 


MATERIAL 

k 

P , 

Kg/a^ 

V/a 

kE^/p 

V-a/kg 

K-Fila 

10 

1750 

2.75x10^ 

4.32x10^“^ 

POLYMIIS 

3.6 

1420 

2.75x10* 

1.9x10^^ 

POLISULPONE 

3.1 

1240 

2.95x10* 

2.17x10^*^ 

Pa.YCARB0NATE 

2.8 

1210 

2.75x10* 

1.75x10^^ 

POLTFROPILENE 

2.2 

900 

8.66x10* 

1.8x10^^ 

TEFLON 

2.1 

2300 

5 .9xl0'^ 

3.2x10^^ 


Substituting values for polypropylene in Equation 58 yields the optiaua 
specific energy of the dielectric which is expressed as: 

w - 7970 (N. 

and for R>20« and (59) 

-0.241 


7970 N' 







The paraneters k and p are material propertiea. For a given dielectric, 
the energy density is directly proportional to the square of the applied 
electric field on the dielectric. To maximize the energy density, the 
dielectric mast be stressed to the mazimnm allowable electric field deter¬ 
mined by the breakdown strength (E^). 

The maximum allowable field is a quantity which depends on the required 
number of cycles, N, (number of charge/discharge cycles), and the voltage 
reversal (R) experienced during the capacitor operation.^ The allowable 
ratio of charging to breakdown voltage reflects the limit on capacitor life 
doe to corona formation.^ When the applied electric field is greater than 
the allowable field, small currents may arise in the dielectric and may 
degrade the dielectric or lead to a complete breakdown.^0 Voltage reversal 
affects the capacitor life because it leads to degradation of the dielectric 
and corona formation. 

Data for high-energy-density capacitors^l were used to determine the 
relationship between capacitor life and the applied electric field. This 
relationship is 


N » (E/E^)-8-3 (56) 

Equation 56 shows that the life of the capacitor is equal to one cycle when 
the applied field equals the breakdown strength and that life increases 
exponentially as the applied field is decreased. 

The effect of voltage reversal on the capacitor life was analyzed in 
Reference 9 and was found to have a negligible effect if the voltage 
reversal is less than 20%. Voltage reversal greater than 20% reduces the 
life by: 


^Sargent, W. J. 'Energy Storage Capacitors', Los Alamos National Scientific 
Laboratory. 

^Dailey, C U, White, C. W., 'Capacitors for Aircraft High Power', AFAH,- 
TR-74-79. 

^^nnhardt, E. E., 'Electric and Magnetic Properties of Materials', Texas 
Tech University. 

^^Naxwell, Data Sheets - Series M Capacitors 
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conditioners will give the possible ES/PC concepts. This is illustrated in 
Table 7. The type of input and output conditioner is also constrained by 
the priaary power source selected and the thruster (e.g.. Series/Parallel 
switching of the primary power sources is only feasible when modular power 
sources are used). No input conditioner may be required when thesmocycles 
such as Rankine or Brayton are the primary power sources and they are used 
with an inertial energy store. The thruster requirements will limit the 
type of energy store and output conditioner which can be selected. 


TABLE 7 

SDMUARI OF ES/PC CONCEPTS 


INPOT 

CONDITIONER 

ENERGY 

STORE 

ODTPUT 

CONDITIONER 

DC/DC Converter 

S/P Switching 

CAPACITOR 

Thyratron 

Spark Gap 

SCR (at closing 
switch) 

LC Circuit 

PFN Circuit 


Normal Inductor 
Superconducting 
Induct or 

SCR (As opening 
switch) 

OC/DC Converter and 
Electric Motor 

0 None 

Coupulsa tor 

SCR (As closing 
switch) 


The next step is to characterize the ES/PC compoikents and develop thi 











INPUT 

CONDITIONER 


Figure 14. ES/PC Concept Using s Nomsl Indnctor as the Energy Store. 


PRIMARY INPUT 

POWER CONDITIONER 


Figure 15. ES/PC Concept Using s Super conduct in 
Store. 


ig Inductors ss the Energ; 


store (Figure 15). A homopolsr generator was also considered but was not 
feasible because of the relatively short pulses required by the thrusters. 


PRIMARY INPUT 

POWER CONDITIONER 


Figure 16. ES/PC (incept Using a Compul sa tor. 


The input conditioners considered for these concepts were: DC/DC 
onverters, series/parallel switching of prisary power modules, and electric 


The ontuut conditioners were switches and inductors. Inducto 


ised concepts (LC and PFN) 
and spark gaps. 
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al application were us« 

to project 

the performance of spa 
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t to power levels 
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i 75 kW. 

The available performance data for 

motors for terre 

str i 1 

appl 

ica tion 

s were curve-fitted t< 

:> obtain the specific power and effi 

cienc 

The 

specif i 

c power is: 






w = 

4.706P0*2926 


(71 

and 

the eff: 

i ciency: 






q = 67,96P0*0252 


(7( 

Ihe 

projected specific power for 

space-qualified moi 

tors is: 




w « 

19.18P®-2368 


(7- 

Thi i 

i Bodel 

hss the ssae ‘slope* 

as the data for 

motors for terre 

strii 


3 .4.4.3 Series-Parallel Switching of Priaary Power Nodoles 
Series-Parallel Switching (S-P) of priaary power aodol 
conditioning method which allows constant power charging of 


be designed to deliwer any required current/vol tage 1 
re.^^ S-P can only be used with aodnlar power sources 
theraoel ectric, and theraionic. 


id Billerbeck, V. J., 'Capacitor Bank Charging by Seri 
ng of Solar Arrays', 16th International Electric Propul 
. AIAA-82-187 9. 











Figure 27. Specific Power of Electric Motor:;. 

To illustrate this couoept, consider a solar array which is used to 
charge a capacitor bank. Ihe solar array can be diwided in N nodules where 
N is an ewen integer. All the nodules are connected in parallel at the 
beginning of the charging cycle (configuration A in Figure 28). The I-V 
characteristic of configuration A is curve 's' in Figure 29. The charging 
proceeds frca point A to point B along carve 'a'. When the voltage reaches 
point B, the nodules are switched to configuration B in Figure 28. The I-V 
characteristic of the array is then curve 'b* in Figure 29. The charging 
proceeds along curve 'b' to point C Siailar switchings are repeated until 
the voltage reaches the required value with the solar array in configuration 
D. Depending on the nunber of nodules and switches, S-P will result in 
capacitor charging which closely approaches the constant power curve, 'e'. 


We assuned that the nodular power source had eight nodules and that thi 
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Figaro 28. Sorieo-Par&llel Switching Concept. 
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Tliis is • total of 21 switches. The mass aad efficiency of each switch is 
calculated separately. The mass of the S-P system is the sum of the 
indiwidnal masses and the efficiency is the product of the individual 
efficiencies. 

3.4.5 Output Conditioner Models 
3.4.3.1 SCR'S 

SCS's are semiconductor switches which have demonstrated long life and 
reliability. The mass, efficiency and operating constraints of SCR's are 
discussed below. 

Fast switching SCR's capable of switching up to 17000 amperes and 2000 
volts in a pulsed mode are available. For higher power requirements, SCR's 
can be placed in a series or parallel arrangement to switch the required 
current and voltage. SCR's have a life higher than the 10^ cycles and an 
efficiency of bettor than 98%. 

SCR current and voltage capabilities are limited by the breakdown 
voltage of the junction and by power dissipation in the on-stste. Assear- 
blies with cooling systems capable of keeping the SCR's at a safe tempera" 
tnre (175® C) are available. 

Manufacturer's specification data^^ were used to model the mass and 
efficiency of the SCR's. Two mass models were generated: one for SCR's 
used as closing switches (pulse rated data were used), and one for S(X's 
used as opening switches (RMS data were used). The mass models are pre¬ 
sented in Figure 30. The mass of SCX's can be best modeled as normalized to 
the mazimiB safe conduction current (m/I). By normal izing the mass, we were 
able to account for the modularity of SCS's. For a specific application, 
the mass of the SCR (and assembly) is the product of the normalized mass and 
the required current. The normalized mass models for the SCS's as closing 
switches are: 

m/I - 1.36 X 10~3 for: V<1200V 

and (78) 

m/I - 0.67 z lO-^V z 5.6 X 10”^ for: V>1200V 


^^Vestinghouse Power Semiconductor User's Manual and Data Book. 
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Figure 30. SCR's Muss Models 
The Bodels for the SCX's as opening switches are: 

m/1 - 2.51 X 10~3 for: V<1200V 

and (79) 

B/I =2.47 X 10"^V -4.5 X 10“^ for: V>1200V 

The efficiency of the SCR's as Bodelled was abont 99%. 

3.4.5.2 Tl^ratrons 

A thyratron is a gas-filled closing switch nsed in aany applications 
which require high peax-power levels (MW or greater) for very short dura¬ 
tions (lO's of ps or less). A thyratron is shown scheBS ti ca 1 1 y in 
Figure 31. It consists of two electrodes: an anode and a cathode. A third 




-ode, the control grid, is placed between the Bain electrodes. These 






Anode 


Sea 1ed 
Envelope 



Control 

Grid 


( 1017 ) 


U-— Cathode 


Figure 31. Tli 7 recron Scliemetic. 

Tabe condaction is controlled by tbe grid voltage. Wlien a positive 
signal is applied to tbe grid, tbe gas in tbe catbode-grid space becomes 
ionized. Ions near tbe grid are then accelerated by tbe anode field, caus¬ 
ing tbe entire volne to become ionized, Tbe tube tben acts as a closed 
svitcb. A tlqrratron will continue to conduct as long as tbe gas remains 
ionized. During conduction, a positive space charge forms around tbe grid 
and effectively removes it from tbe circuit. Additional potentials applied 
to tbe grid will not control tbe anode-to-catbode current. In order to 
restore tbe tube to its 'open' condition, tbe anode voltage must be removed 
or made slightly negative for a time sufficient to allow tbe gas to 
deioniz e. 

Tbyratrons are used with mercury, rare gases, or hydrogen as tbe filler 
gas. Hydrogen is tbe most commonly used filler gas. Hydrogen provides a 
lower recombination (recovery) time and tbe voltage drop across its plasma 
is low enough to keep tbe ion velocities low. Ions with high velocities can 
damage tbe electrodes. 

Tbe paraieters affecting tbyratron performance are: peak current, peak 
voltage, pulse width, recovery time, pulse frequency (rep-rate), forward 
voltage drop, and lifetime (total niaber of pulses), 

Tbe peak current is limited by cathode arcing and grid quenching. When 
tbe peak current exceeds some allowable limit, an arc forms between tbe 
anode or grid and tbe cathode and damages the cathode. High peak currents 
can also cause localized beating (only a small part of tbe cathode is used 
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for conduction in each pulse) which degrades the perforaiance of the thyra- 
tron. Thyratrons with a peak current capacity of 120 kA have been buil t.^^ 

The peak voltage which a thyratron is capable of 'standing off is 
detemined by the distance between the anode and the control grid. Under 
repetitive pulse operation, more complex mechanisms begin to limit the 
stand-off voltage, such as charge accumulation on the tube walls. One way 
to increase the stand-off capabilities is to increase the number of control 
grids. Manufacturing procedures and tube packaging constraints place a 
practical limit of 40-50 kV/gap on the stand-off voltage. 

The pulse width is limited by the peak current and the desired life of 
the tube. High currents for long times can result in arc formation and a 
high rate of cathode vaporization. Typical pulse lengths at peak currents 
are on the order of tens of microseconds. 

The pulse rep-rate and the recovery time are inter—relaved. The maxi¬ 
mum rep-rate is limited by the recovery time (time required for a gas to 
recombine and regain its stand-off capabilities). Thyratrons with rep-rates 
on the order of kHz have been built. High rep-rates increase the average 
heating of the tube and affect the tube performance and life. 

The voltage drop in hydrogen thyratrons must be limited to keep the 
ion velocity below the critical velocity which damages the electrodes. 
Typical voltage drops are about 100 volts. 

Thyratron life is strongly related to the allowable peak current, pulse 
width, and rep-rate in high frequency applications. The end of the tube 
life occurs when the cathode can no longer emit electrons into the neutral 
plasma to begin conduction. Typical hydrogen thyratron life is greater than 
lO-^l pulses. 

A mass model of the hydrogen thyratron was developed by curve-fitting 
available data. The data was derived from manufacturers' specification 
sheets. The mass of thyratrons can best be presented as specific power (w) 
versus peak power (P). The relationship is: 

w = 5.31 X 108 p-0-21 (go) 


and G Ceramic Metal Grounded Grid Thyratrons Data Sheet, H5003-1. 
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two electrodes (an anode and cathode) separated by some distance (called the 
inter-electrode spacing). The elements are enclosed in a sealed envelope 
which is either evacuated or filled with a pressurized gas. 

The direct over—vol ted spark gap breaks down when the stress in the 
inter-electrode gap exceeds some critical value. This critical value is 
called the self-breakdown voltage (SBV). The self-breakdown voltage is a 
very strong function of the inter-electrode spacing, the electrode shape, 
and the filler gas conditions (if a gas is present). 

The triggered spark gap has a third electrode placed inside, but sepa¬ 
rate from, one of the main electrodes. An electrical pulse is applied to 
the third electrode generating an arc between it and the nearby electrode. 
The resultant 'mini discharge' causes the main electrodes to break down. 
The breakdown time is a function of the voltage across the gap. The closer 
the gap voltage is to the SBV, the faster the gap will break down. 

The field distortion gap has a third element (usually a ring shaped 
device) placed between the two main electrodes. This 'trigger' is placed at 
a point where it is not normally distorting the electric field lines in the 
gap. When an electrical pulse is applied to the trigger electrode, the 
field in the gap is distorted and a breakdown occurs. 

The parameters which affect spark gap performance are: peak current, 
peak voltage, pulse width, pulse rep-rate, total coulomb transfer, recovery 
time, and lifetime.^ 

The peak current is limited by the required pulse width, the resultant 
charge transfer per pulse, and the total niaber of pulses. Peak currents of 
hundreds of kA are not uncommon in typical spark gap applications. 

The voltage which a spark gap can stand off is a complex function of 
the filler gas type (if a gas is present), the electrode material, and the 
electrode spacing. For extremely high voltage applications, the siaximimi 
standoff voltage is, in practice, limited by surface flashover outside of 
the gap rather than by the gap design itself. Peak voltages of millions of 
volts have been achieved and values as high as 100 kV are typical. 

The maxim» pulse width that a spark gap is capable of is a function of 
the required current and life. Pulse widths of tens of microseconds are 
common. The maximum single-shot pulse width is limited by electrode 
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cooling. The minimma pulse width is sot by external circuit parameters. 
Pulse widths less than 1 ns have been achiewed. 

The allowable pulse rep-rate is limited by the recovery time and the 
rated voltage standoff. The recovery time is a fundamental limit to the 
maximoB pulse rep—rate. The recovery time is limited by the plasma recombi¬ 
nation time and the gas and electrode cooling rates. For 100% voltage 
standoff, a practical limit on pulse rep-rate is a few hundred pulses per 
se cond. 

The total charge transferred during the life of a switch is equal to 
the charge transferred per pulse multiplied by the total nianber of pulses. 
The charge transferred per pulse is directly related to the current and 
pulse width. Spark gaps transfer a relatively large charge per pulse. This 
changes the surface microstructure of the main electrodes and affects 
operation. Very large charges per pulse have been transferred in single¬ 
shot devices. 

The life of spark gaps is an area of considerable uncertainty. Nearly 
all manufactiirers estimate spark gap life at 10,000 pulses at maximum rated 
conditions. Very little work on extrapolation of life data has been per¬ 
formed. The charge transferred per pulse must be decreased to increase life 
(a spark gap can only transfer a finite total charge during its life). 
Increased life can also be obtained by preionizing the gas with an external 
source to reduce the charge transfer load imposed on the electrodes. Al¬ 
though this technique has not been thoroughly tested, it is expected to 
increase the life to greater than 10^ pulses. 

Many factors affect the performance and mass of spark gaps. After 
analyzing the available data, we determined that the mass of the spark gaps 
is primarily a function of the total charge transferred. With this 
approach, we account for the peak current, pulse width, and life of the 
spark gap. The resulting mass model is presented in Figure 33 and is 
written as: 


m = 0.0525 qOAAII 


(83 ) 


where m is the mass in kg and C is the total charge transferred in coulombs. 
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4.2,3 Veloc: 


ity Increment 

Uie optiaiBD propulsion systen mass was proportional to the vel 
increment for s ta t' onkeeping missions (Figure 41). We assumed tha 
velocity increment for orbit transfer missions is constant. 
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Figire 41. Velocity Increment Effects on Optimum Propulsion System Ma 

4.2.4 Primary Power Options 

Photovoltaic power sources resulted in the lowest mass of og 
propulsion system for stationkeeping systems. This is illustrat 
Figure 42. 

Nuclear-based primary power sources become competitive with 
voltaic for orbit transfer missions. A typical example is illustrai 
Figure 43, 

4.2.5 ES/PC Options 

The capacitive-based ES/PC options resulted in the lowest-mass o 
propulsion systems for stationkeeping missions. The normal inductive 
ES/PC options resulted in comparable optimn propulsion system mass. 
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la Table 12 we are preseating two propnlsion systens designed by 
EPSDES. The first has an Isp of 1600 s. similar to the Martin Marietta 
system, and the second has an Isp of 3100 s. The propnlsion system designed 
by EPSDES to operate at an Isp of 1600 s has the same mass as the Martin 
Marietta system. We cannot compare other characteristics of the system 
becanse of lack of information on Martin Marietta's system. But, we believe 
that the EPSI£S-se 1 ected system operates with lower pulse mass and higher 
frequency than Martin Marietta's. The second system presented in Table 12 
is the optimimi system selected by EKISS. The optimum system has 30% less 
mass than the Martin Marietta system and operates at an Isp of 3100 s. We 
believe that the optimnn system has a higher frequency and lower pulse mass 
than Martin Marietta's system. 

4.2 SOISITIVITT ANALYSIS 

We analyzed the sensitivity of optimmi propnlsion system mass to varia¬ 
tions in mission parameters, primary power options and ES/PC options for 
stationkeeping and orbit transfer missions. This was done by exercising 
EPSDES over the entire range of interest. 

4.2.1 Payload Mass 

The optimum propulsion system mass was proportional to payload mass for 
both stationkeeping and orbit transfer missions. This is illustrated in 
Figures 37 and 38 respectively. The mass of the propnlsion systems for 
stationkeeping missions ranged from 5 to 15% of the payload mass depending 
on the thruster used. The EBG had the highest propulsion system mass and 
the other thrusters had similar masses. Ihe mass of the propnlsion systems 
for orbit transfer missions ranged from 2.5 to 3.5 timet the payload mass. 
Ihe EBG had the highest propulsion system mass. 

4.2.2 Mission Duration 

The optimum propulsion system mass was proportional to mission duration 
for stationkeeping missions (Figure 39), and increased exponentially with 
decreasing mission duration for orbit transfer missions (Figure 40). The 
optimum exhaust velocity decreased with decreasing mission durations for 
orbit transfer missions. 
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TABLE 11 

ORBIT mANSFER PERFORMANCE COMPARISON 


















TABLE 10 


STAIIONKEEPING PERFORMANCE COMPARISON 
(Hie specific energy of cepacitors was assmed to be 1.7 J/kg) 


Paraaeter 

lAP MPD 
Systea 
N-S 
Module 

lAP MPD 
Systea 
E-W 
Module 

lAP MPD 
Systea 

MM^S MPD 
Systea 

Chosical 

Systea 







Payload Mass, kg 

20978. 

20978. 

19878 

19878 

19 878 

AV, a/s/yr 

30.5 

19.5 

100. 

100. 

100. 

Mission Tiae, yr 

5. 

5. 

5. 

5. 

5. 

Thruster No. 

1 

1 

4 

19 


Average Power, kw 

6.39 

7.2 

27.2 

4.9 


Thruster Eff., % 

11.28 

11.35 

11.3 

54. 


Total Eff. . % 

9.79 

9.73 

9.76 

43. 


Isp. s 

4000. 

3500. 

3700. 

5000. 

285. 

Frequency, hr 

80.6 

73.7 

77.1 

2.7 


Pulse Energy, J 

69. 

83.7 


625. 


Pulse Tiae, as 

0.1 

0.11 


1. 


ES/PC Mass, kg 

102. 

114. 

432. 

1476. 


Primary Power, kg 

121. 

132. 

506. 



Dry Mass, kg 

239. 

258. 

994. 

2579. 

1142. 

Propellant Mass,kg 

81. 

59. 

280. 

233. 

3984. 

Propulsion Mass,kg 

320. 

317. 

1274 

2812. 

5126. 


Hie EPSOES selected a propulsion systea that has a mass less than 1/2 
the Martin Marietta systea. Ihe selected systea has an Isp of 3700 s which 
is twice the Isp of the Martin Marietta systea. This led to higher systea 
efficiency and lower aass. 

We suspect that the higher Isp was selected because of a difference in 
the specific energy of the energy store (capacitor). Therefore, we reran 
EPSDES with a specific energy of 1.7 J/kg which we believe that Martin 
Marietta used in their analysis. The results of the new run of EPSDES and 


Martin Marietta results a: 


presented in Table 12. 
















4.1.3 .1 Hie Stetionkeepi&g System 

The stationkeeping mission specifications were: 

Payload mass - 19878 kg 

Mission time = 5 years 

AV = 500 m/s 

Martin Marietta assmed fonr thrusting modules. Three modules had five 
thrusters and one module had four thrusters. Each module had one ES/PC 
systeiB. N-S stationkeeping required a AV of 61 m/s/yr and a duty cycle of 
0.666. E-W stationkeeping required a AV of 39 m/s/yr and a duty cycle of 
0.333. For EPSDES modeling purposes, we assimied that two modules were used 
for N-S and two for E-V stationkeeping, that each module had one thruster, 
and that each thruster was required to deliver half the required AV. The 
payload mass was also increased by a mass equivalent to 3 thrusting modules. 
Using EPSDES computer code, both E-W and N-S MPD propulsion systems were 
sized. The EPSDES and Martin Marietta results are presented in Table 9. 

The system selected by the EPSDES has a propulsion system mass of less 
than 1/3 of the Martin Marietta system. The main difference between EPSDES 
and the Martin Marietta results was the mass of the ES/PC EPS1£S predicted 
an ES/PC which was one-fifth the mass of the ES/PC predicted by Martin 
Marietta. This was caused by the conservative specific energy of capacitors 
used by Martin Ma v'tta, which was 1.7 J/kg. As a result of this compari¬ 
son. we reran EPSDES for the same mission but we assmed that the capacitors 
had the same specific energy used by Martin Marietta. The EPSDES results of 
the new run and Martin Marietta results are presented in Table 10. 

The new system selected by the EPSDES has a propulsion system mass of 
less than a half of the Martin Marietta system. EPSEES selected a propul¬ 
sion system that delivers small pulses one-eighth the size of Martin 
Marietta's at a lot nigher frequency than Martin Marietta (25 times). This 








4.1.2 OptiBnm Propalsion Systems 


We used EPSOES to examine tke general characteristics of optimized 
propalsion systems for both stationkeeping and orbit transfer. The mass of 
the optimm propalsion system for stationkeeping missions ranges from 5% to 
15% of the payload mass depending on the thrnster ased. Ihe EBG had the 
highest propalsion system mass. The minimnm al lovable pal se mass of the EBG 
is 0.1 g, 3 orders of magnitade higher than that for the other thrasters. 
This resnlts in high stored energy and a corresponding high ES/PC mass. The 
largest mass component of the propalsion system for sta tionkeeping was the 
propellant at about 50%. The optimization tradeoff for stationkeeping is 
primarily between energy storage mass in the ES/PC and primary power mass. 
As expected the ES/PC and primary power have about equal mass in optimized 
systems. 

The mass of the optimm propalsion system for orbit transfer missions 
ranges from 2.5 to 3.5 times the payload mass, with the EBG being the 
heaviest. The ES/PC mass was the least massive component at about 15% of 
propalsion system mass. In orbit transfer missions the optimization 
tradeoff is primarily between propellant and primary power. As expected 
these two components had about equal mass in optimized orbit transfer 
systems. 

Within the ES/PC system, the major components of optimized station¬ 
keeping systems had approximately equal mass. The EBG was an exception. The 
major mass component of the ES/PC for the EBG was the energy store at about 
80% of the ES/PC mass. The major mass component of the ES/PC for orbit 
transfer missions was the radiator at about 45%. The energy store was 30%, 
the input conditioner 1S%, and the output conditioner 10% of ES/PC mass. 

4.1.3 Comparison of EPSDES with Martin Marietta Results 

Martin Marietta^^ conducted two propulsion system point designs in 
their 'MPD Definition Study Program'. The first point design was for a 
stationkeeping mission and the second was an orbit transfer mission. We 
exercised EPSDES for the same missions. Both designs were MPD thruster- 
based. Our objective was to compare the EPSDES approach to system sizing 
with the more traditional approach used by Martin Marietta. 

^*'MPD Thruster Definition Study,' Martin Marietta, 1983. 
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Figure 35. Propulsion Systea Optimization (Orbit Transfer, TPP). 
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SECTION 4 


RESULTS 

We exercised the EPSDES progrta extensively. The first series of 
conputstions were designed to exaaine the perfozBsnce of the progrsa, to 
establish the validity of the solutions which it finds, and to detemine the 
general characteristics of optinally sized propulsion systems. We then used 
EPSI£S to conduct a sensitivity analysis. We examined the sensitivity of 
propulsion system size to variations in mission parameters over the entire 
range of mission parameters defined in Section 2. 

4.1 EPSDES VALIDATION 

4.1.1 EPSDES Optimization 

The basic assomption underlying the development of EPSDES was that for 
each mission and thruster combination, there is one set of propulsion system 
output parameters (Up, and f) which minimizes the propulsion system mass. 
The first task was to verify this asstBption. 

We selected an orbit transfer mission and a TPP-based propulsion system 
to check the validity of EPSDES. EPSDES was run with a fixed exhaust 
velocity and the pulse mass (and frequency) was varied. The results at 
several exhaust velocities are shown in Figure 35. Each point on the graph 
represents an integrated TPP-based propulsion system which will be capable 
of doing the required mission. There is a clear minimum in propulsion 
system mass at each exhaust velocity. We conclude that for each exhaust 
velocity, there is an optima pulse mass which minimizes the propulsion 
system mass. 

Figure 35 also shows that selection of the exhaust velocity has a 
pronounced effect on propulsion system mass, ^e locus of the minima of the 
curves in Figure 35 are plotted in Figure 36. There is a combination of 
exhaust velocity and pulse mass (and frequency) which results in minimum 
propulsion system mass. This combination results in the 'optimized' propul¬ 
sion system. We ran EPSDES for a wide range of missions and an optimum 
propulsion system was always found. 
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3.5 SYSTEM INTEGHATION 


To assess the propalsion system as an integrated system, the subsystems 
and components characteristics and performance models must be incorporated 
into a system level description of an electromagnetic propulsion system. 
The system integration can be best accomplished by using numerical 
techniques because of the complexity of the models and the large ntimber of 
options. 

We developed an Electromagnetic Propulsion System Des ign (EPSDES) com¬ 
puter program to assist in selecting and sizing the optimum propulsion 
system for a specific mission. EPSIKS is written in FORIXAN. EPS1£S calcu¬ 
lates the combination of m^, Vq< ^nd f which minimizes the propulsion system 
mass and satisfies the mission requirements. A simplified flow chart is 
presented in Figure 34. 

EPSDES is an interactive program which prompts the user for the re¬ 
quired input data. The input data are: (1) the mission parameters AV, 

tj,, and a, (2) the thruster type, (3) the primary power type, (4) the ES/PC 
type, (5) the exhaust velocity range and increment, and (6) an initial guess 
for the pulse mass. The output data are the optirntmi subsystem masses and 
component characteristics at each exhaust velocity within the range 
specif ied. 

The validation of EPSDES as a propulsion system optimization program is 
presented in Section 4. Also, we are presenting in Section 4 a comparison 
between EPSDES results and known point designs, and the results of sensi¬ 
tivity analyses of the optimum propulsion system size to variations in 
mission parameters. 
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superconductor inductive ES/PC resulted in a much heavier optimuin propulsion 
system. The inertial encrev stores were not feasible for the PIT and TPP 



Exhaust Velocity (Um/s) 

Figure 44. Propulsion System Mass for Different ES/PC Options for 
Stationkeeping Missions. 

For orbit transfer missions, the capacitive-based ES/PC options re¬ 
sulted in the least massive optimum propulsion system. The inductive-based 
ES/PC options (both normal and superconductor) resulted in comparable re¬ 
sults. The inertial energy stores resulted in much more massive propulsion 
systems. The propulsion system masses for various ES/PC options are pre¬ 
sented in Figure 45. 


4.3 POINT DESIGNS 

To demonstrate the methodology developed in this program, we conducted 
two point designs, one was for a stationkeeping mission and one for an orbit 
transfer mission. The missions specifications and the results of the point 
desiens are discussed below. 
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propulsion system is presented in Table 13. the size and cbaracteristics of 
tbe thrusters in Table 14, and the size and characteristics of the ES/PC 
components in Table 15. 

The optimtB propulsion system output parsmeters are: a pulse mass of 
OJ mg, an ezhaust velocity of 50,000 m/s. and a frequency of 7 Hz. These 
output pameters resulted in a propulsion system which has a mass of 1200 
kg, which is 6% of the payload mass, and a system efficiency of 22%. The 
average power level required to power the system is 19 kW. 

The low system efficiency is caused by the low efficiency of the 
thrusters which is about 27%. The thrusters have a low efficiency because 
of the low pulse duration which is on the order of few hundred microseconds. 

The optimal ES/PC mass is 350 kg which is 30% of the propulsion system 
mass. As expected, the major mass components of the ES/PC have equal 
weights. 


TABLE 13 

STATIONBEEPING OPTIIOM HIOPHLSION SYSTEM FOR A POINT DESIGN 


Parameter 

MH) 

System 

N-S 

Module 

MPD 

System 

E-W 

Module 

Propulsion 

System 

Payload Mass, kg 

21000.* 

21000.* 

20000. 

A V, m/s/yr 

31.25 

31.25 

125. 

Mission Time, yr 

7. 

7. 

7. 

Thrust/Engine, N 

0.07 

0.021 

0.045 

Average Power, kw 

6.27 

2.80 

18.14 

Ihrust Eff. , % 

30.90 

23.06 

26.98 

Total Eff. . % 

25.52 

19.00 

22.25 

Isp, s 

4500. 

5000. 

4750. 

Frequency, Hz 

5.40 

6.97 

6.18 

Pul se Energy, J 

958. 

331. 

644.5 

Pulse Time, ms 

0.31 

0.18 

0.24 

ES/PC Mass, kg 

110. 

62.7 

345.4 

Dry Biass, kg 

250.7 

157.5 

816.4 

Propellant Mass,kg 

103.3 

92.5 

391.6 

Propulsion Mass,kg 

354. 

250. 

1208. 


Ihe payload mass was increased by an amount equivalent to the mass of 3 
thrusting modules. 
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TABLE 14 


ISBDSTER SIZE AND CHAEACTERISTICS FOR STATIONEEEPING POINT I^SIGN 


Parameter 

MPD 

for 

N-S 

MPD 

for 

E-f 

Outer radius, cm 

2.8 

1.9 

Length, cm 

11.0 

5.9 

Efficiency. % 

30.9 

23. 

Mass flow rate, g/s 

0.95 

0.33 

Isp. s 

4500. 

5000. 

Frequency. Hz 

5.4 

6.97 

Pulse duration, ms 

0.31 

0.18 

Pulse energy. J 

957. 

330.7 

Average Voltage. V 

196.5 

183. 

Average current. A 

15810. 

9889. 


4.3.2 Orbit Transfer Misaioa 

Two thrusters were sssuned for this mission. The required AV was 
divided equally among the two. 

The optimm propulsion system for the orbit transfer mission is pre¬ 
sented in Table 16. the size and characteristics of the thruster in 
Table 17. and the ES/PC system in Table 18. 

The selected ES/PC system was a DC-DC converter, a capacitor and a 
hydrogen thyratron closing switch. 

The optimna propulsion system output parameters are: a pulse mass of 
35 mg. an exhaust velocity of 26000 m/s and a frequency of 10.5 Hz. The 
propulsion system mass is 20.400 kg which is twice the payload mass. The 
system efficiency is 55%. Ihe average power level required to power this 
mission is 450 kW. 

The mass breakdown of the major components of the propulsion system is 
as follows: the primary power source mass is 10.000 kg. the propellant and 
containers mass are 7660 kg. and the ES/PC mass is 2600 kg. The primary 
power mass is the heaviest because of the relatively short mission time. 

The heaviest component of the ES/PC was the radiator with a mass of 
865 kg (35% of the ES/PC mass). Ihe energy store mass was about 30% of the 
ES/PC mass, and the input and output conditioners have about equal mass (15% 
of the ES/PC mass). 
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TABLE 16 


ORBIT TRANSFER PROPULSION SYSTEM 


Parameter 

Propulsion 
System 
for 1 PIT 

Total 

Propnlsi on 
System 

Payload Mass, kg 

20000.* 

10000. 

AV, m/s 

2900. 

5880. 

Mission Time, days 

100. 

100. 

Thrust/Engine, N 

9.60 

9.60 

Average Power, Mf 

225. 

450. 

Thrust Eff.. * 

62.3 

62.3 

Total Eff. , % 

55.5 

55.5 

Isp, s 

2600. 

2600. 

Frequency, Hz 

10.5 

10.5 

Pulse Energy. kJ 

19. 

19. 

Pulse Time, ps 

16.5 

16.5 

ES/PC Mass, kg 

1300. 

2600. 

Dry Mass, kg 

7010. 

14020. 

Propellant Mass,kg 

3190. 

63 80. 

Propulsion Mass.kg 

10200. 

20400. 


The payload mass was increased by an amount equivalent to the mass of 
one thrusting module. 


TABLE 17 

THRUSTER SIZE AND CHARACTERISTICS FC« ORBIT TRANSFER MISSION 


Parameter 

PIT 

Diameter, m 

2.18 

Efficiency, % 

62.3 

Isp, $ 

2600. 

Frequency, hz 

10.5 

Pulse duration, ps 

16.5 

Pulse energy, kJ 

19. 

Average Voltage, kV 

24.6 

Average current, kA 

126.6 
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TABLE 18 


ES/PC SIZE AND CHARACTERISTICS FOR ORBIT TRANSFER MISSION 


Compoae&t 


Capacitor 


Paraae ter 


Efficiency, % 

Mass, kg 
Volnme, m3 
Dimension, m 
Capacitance, pF 
Voltage, kV 
Energy stored, kJ 
Number of cycles 
Specific energy, J/kg 


[ ES/PC 
System 
for 1 PIT 

98. 

402.4 

0.22 

|0.61 

62.7 

124.6 

ll9. 

[9.1E7 

47.3 


Tbyra tron 


Efficiency, % 

Mass, kg 
Peak power, W 
Pulse Duration, ps 
Current, kA 


98.9 

174.8 

1.16E9 

16.5 

126.6 


DC-DC 
Converter 


Efficiency, % 
Mass, kg 

Average power, kW 
Volume, m3 
Dimension, m 


92. 

t204. 

1225. 

I7.5E-3 

|0.42 


Radiator 


Temperature, Deg. E 
Mass, kg 

Surface area, m^ 

Avg. power rejected. 


kw| 


300. 

432.5 

134.1 

24.3 
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SECTION 5 


CONCLUSIONS AND RECOMMENDATIONS 


S.l CONCLUSIONS 

A nnaber of specific conclnsions cen be drnrn from tbis stad^. 

(1) Optiaizstioa of integrated propulsion systesis requires that the prinary 
pover sources, the ES/PC, and the thruster perfomance be related to 
the primary propulsion system outputs (pulse mass, exhaust velocity, 
frequency, and total thrusting time). 

(2) Models which meet these requirements were successfully developed for 
all the systems of interest. 

(3) Most of the models developed were based on an understanding of the 
physical principles of operation and a nirntber of assumptions concerning 
the performance of optimized systems. The rest of the models were 
developed by curve-fitting available perfomance data, either from 
manufacturer's literature or from point design studies. 

(4) The excellent qualitative and quantitative agreements between model 
prediction and observed performance lend considerable credence to the 
models. 

(5) The models are extremely versatile, pemitting rapid evaluation of 
proposed alternatives fox specific missions. 

(6) The models can easily be updated when new data becomes available. 

Using these models to analyze optimized propulsion systems for sta- 

tionheeping and orbit transfer missions, we can conclude: 

(1) There exists one set of propulsion system outputs that minimizes the 
mass of the system. 

(2) The optimum propulsion system outputs are within the capabilities of 
the state-of-the-art thrusters. 

(3) The propulsion system mass is about 10 percent of the payload for 
stationheeping and ranges from 2.5 to 3.5 times the payload mass for 
orbit transfer missions. 

(4) Photovoltaic power sources are the least massive for stationkeeping 
missions and nuclear-based power sources become competitive for orbit 
transfer missions. 

(5) Average power levels range from 2 to 35 kW for stationiceeping and from 
200 kW to 2 MW for orbit transfer missions. 

(6) Capacitive based ES/PC systems are the least massive, inductive based 
systems are comparable, and inertial systems are much heavier. 
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(7) ES/PC conponents have equal weights for stationheeping systems, and the 

radiator is the largest component in the system for orbit transfer 

missions (45%). 

5.2 RECOMMQIOAXIONS 

The thruster performance plays an important role in sizing the ES/PC 
and the primary power source. It is necessary to validate the performance 
models of the thrusters over the wide range m^ ^ad v^ of interest. There 
are two major areas of uncertainty in the thruster models developed. 'Die 
first problem arises from the nature of the thrusters for which data was 
available. It is not clear that the experimental thruster perfomance data 
available to us was generated with optimimi thrusters (that is, the thruster 
was optimized for every and v^ for which data was reported). The per¬ 
formance data must be evaluated more carefully to answer this question. The 
second problem arises from the limited amount of data available. The 
missions considered in this study require the thruster models apply at 
operating conditions (m^ and v^) which could be orders different from the 
conditions for which we have data. Therefore, we recommend that performance 
data for optimum thrusters be generated over a much more extensive range of 
Bp (or mdotp) and v^. 

The effects of the component life (nmiber of cycles) on most of the 
components forming the ES/PC are not well understood. The perfomance of 
the components was degraded to account for the component life effects. We 
recommend that a more extensive analysis be conducted to evaluate the life 
effects on the components performance of the most promising ES/PC. This can 
be done by first selecting the most promising ES/PC Our analysis shoe that 
the capacitive—based ES/PC is the most promising. Then each component in 
the ES/PC should be tested individually to study the life effects on the 
performance of the component. After testing the components individually, 
they should be integrated and tested as a system to insure that no com¬ 
pounded effect is available. 
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APPENDIX A 


PRIMARY POWER FOR ELECIROMAGNETIC PROPULSION 

Priaftry power sonroe perfoxaence and aata aodels are required to 
select, size, and design an integrated electroaagnetic propulsion systea. A 
priaary power source delivers electrical output at a relatiwely low average 
level to an Energy Storage/Power Conditioning systea (ES/PC) vRicli condi¬ 
tions the power and delivers high-power pulses to the electrical thruster. 
The priaary power source is defined as the 'raw' power sonroe and all the 
coaponents (radiator, generator, converter, etc.) required to deliver 
electrical power to the ES/PC 

The perforaance characteristics and aass aodels for the priaary power 
sources were detezained by conducting a literature survey to gather the 
necessary data, contacting technical experts to confim the results and 
asking refineaents where needed. The literature contains considerable point 
design data but very little pametric or 1 iait data. 

A wide range of priaary power options were considered including: solar 
cell arrays, solar theraophotovoltalc, nuclear theraionic, nuclear thezao- 
electric, nuclear Rankine cycle, nuclear Brsyton cycle, and nuclear MED. 
Pariaetric perforaance aodels were developed for power levels ranging froa a 
few kilowatts to 10 MW. 

A.1 SOLAR (£LL ARRAYS 

Solar cell arrays are designed in collapsible sections and are opti- 
aized to ainiaize power transaission losses. Transaission losses are olaic 
losses and are ainiaized by power transaission at the highest possible 
voltage. Olaic losses are directly proportional to the current squared for 
a fixed bus bar cross-sectional area.^ 

The perforaance of the solar cell arrays is affected by: (1) plasaa 
arcing, (2) voltage breakdown, (3) voltage transients and Van Allen Belts. 


^Baner, David P., Barber, John P., Swift, Hal lock F.. Vahlberg, C Julian, 
'Electric Rail Gun Propulsion Study (Advanced Electric Technology—High 
Thrust).' AFRPL-TR-81-02, 1981. 
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Solar cell arrays degrade when passing through the Van Allen belts. Degra¬ 
dation of up to 40% is reported.^ Provisions anst be aade for this degra¬ 
dation if solar cell arrays are to be nsed in orbit transfer missions. 

Solar cell arrays can be designed to deliver ary power level required 
since they are composed of solar cells oonnected in series and parallel. 
Solar cell arrays capable of providing power levels up to 10 kW have been 
deployed in space (nsed on Skylab). Arrays capable of multi-hundred kilo¬ 
watt have been designed and analyzed. 

Operating voltages ranging from 28 V to 900 V have been reported. 
Keferenoes 4. 5. 6 , 7, and 8 report the performance of arrays delivering 
multi-hundred watts with operating voltage ranging from 28 to 42 VDC. Ref¬ 
erences 7 and 9 report analyses of solar power systems ranging from 100 to 
300 V. References 10 and 11 reported tradeoff analyses for power management 
and transmission of multi-hundred kilowatt and multi-megawatt systems. Hie 
multi-hundred kilowatt system had an operating voltage of 750 for the DC 
transmission and an operating voltage of 440 for an AC/DC transmission. 


^Taussig. R., et al., 'Overview Study of Space Power Technologies for the 
Advanced Energetics Program.' NASA CR-ld5269. 1981. 

^Brooks, G. R., 'Orbiting Solar Observatory (OSO-8) Solar Panel Design and 
in Orbit Performance.' Intersociety Energy Conversion Engineering Conference 
(lECEC). 789321. 1978. 

^Thronton. 'AIM Solar Amy In-Flight Performance Analysis.' lECEC. 1974. 

^Goldhammer. L. J.. 'Design and Flight Performance of the Pioneer Venus 
Multiprobe and Orbiter Solar Arrsys.'IECEC. 809076. 1980. 

^Lnkens. F. E.. 'Advanced Development of a Programmable Power Processor,' 
lECEC. 809076. 1980. 

^McKinney, H. N., Briggs. 0. C., 'Electrical Power Subsystem for the 
Intelsat V Satellite.' lECBC. 789084, 1978. 

^Glass, M. C., 'A Six Kilowatt Transformer-Coupled Converter for Space 
Shuttle Solar Power Systems.' lECEC. 809146. 1980. 

^®Mildice, J. W., 'Study of Power Management Technology for Orbital Mul ti- 
100 kf Applications,' NASA CR1S9834. July 1980. 

^^Peterson, D. M. and Pleasant, R. L., 'Study of Mul timegawatt Technology 
Needs for Photovoltaic Space Power Systems,' GDC-AST-81-019, August 1981. 







lannciies tad coastraction ia space. 


System mass data gatliered from refereaces 1-18 are preseated ia 
Figure A-2. Figure A-2 shovs tliat the specific mass is about SO kg/kW for 
oae kilowatt arrays sad decreases asymptotically to 10 kg/kW for 100 kW sad 
above. This data was curve-fitted to geaerate a mass model: 

■,p « 0.013 + 40.61/Ppp (A-1) 

where m^^ ia the specific mass ia kg/W aad Ppp i* the average power level ia 
watts. Equatioa A-1 is true for power levels above 1000 W. For power 
levels lover thaa 1000 W. the mass of the system is assumed to be 54 kg 
iadepeudeat of power level. Aa overall degradatioa of 20% ia the perfor- 
maace of the arrays was assumed for orbit traasfer missioa. 

A.2 SOLAR IHERJfOF&OTOVOLTAIC 

The thermophotovol taic (TPV) coacept offers a viable alternative to 
coBveational solar cell arrays. It offers a high efficieacy (50%),thus 

^^Briggs, 0. C, Pollard, H. E., 'Intelsat V Solar Array Design aad Develop¬ 
ment Summary,' lEdC, 7 8 9322, 1978. 

^^Kelly, F. G., Luf t, W., Kurland, R. M., 'Design Features of the TEDSS 
Solar Array,' lECEC, 78 9323, 1978. 

^^Borduas, H. F., et al., 'Medium Power Deployable Hybrid Solar Arrays,' 
lECEC, 780328, 1978. 

^^Peterson, 0. G., 'INSAT-1 Solar Array Design aad Development Summary,' 
lECEC, 19 80. 

^^Souza, C J., 'Large Area Flexible Solar Array Design for Space Shuttle 
Application,' lECEC, 809485, 1980. 

^^Woodcock, W. G. Ill, 'Multi-Hundred Kilowatt Solar Arrays for Space,' 
lECEC, 809271, 19 80. 

^^Trumble. T. M., 'Space Application for Gallium Arsenide Solar Cells,' 
lECEC, 82 9 269, 1982. 
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Figare A-2. Specific Mess of Soler Cell Array Systems. 

a redaction in photoTOltaic cell area and a possible redaction in weight. 
The photowoltaic cells are shielded from the external enwirossoent and do not 
soffer radiation degradation to the extent that conventional panels do. 

TPV power systems can be designed to deliver any power level required 
and the same limiting factors which apply to solar cell arrays (voltage 
levels, breakdown, and transients) apply to TPV systems. 

IPV systems are still in their infancy and many parameters which affect 
performance are still not well onderstood. Design criteria and methods have 
not been developed. A literature survey revealed that feasibility and 
sensitivity studies have been conducted.^®Very little specific mass 
data was located. Reference 19 reported some specific mats data, which we 
used in conjunction with a system energy balance model (Figare A-3), 


^^Swanson. R. M., et al.. 'Silicon Photovoltaic Cells 
Conversion,' EPI-ER-478 Progress Report, Stanford Ele( 
1977. 


1 Thermophotovol taic 
onics Lab., February 


'‘^Horne, V. E., et al., 'Improved Thermophotovol taic Power Systems,' lECEC, 
82 9018, 1982. 



















The assomptions made in order to model the specific mass for large 
arrays were: 

(1) 72% of the incident solar energy is retained in the TPV converter. 

(2) Ihe efficiency of the photovoltaic cells is 40%. 

(3) One concentrator is nsed if the required power level can be 
achieved with a concentrator radins less than 3m. (limitations of 
the available data), and mnlti-3m. radins concentrators are used 
for higher power levels. i.e., the specific mass is constant for 
power levels requiring more than 3m. radins concentrators. 

We assumed that large TPV panels will be formed with large concen¬ 
trators, eventhongh smaller concentrators have a lower specific 
mass, because we felt that their advantages outweigh the possible 
savings in mass of the TPV system. The advantages of the large 
concentrators are: they minimize the number of concentrators 
required to achieve a power level (leads to less complex system), 
they reduce mantifacturing slope errors of concentrator surfaces, 
and they simplify the sun-TPV alignment requirement. 

The specific mass model for power levels ranges from 1 hW to 11 kW is: 

»sp * 0.005 + 4.6xlO"‘^Ppp (A-2) 

where m^^ is in kg/W and Ppp i* the average power level in watts. For 
higher power levels the specific mass was assumed constant and equal to lOJ. 
kg/kw. Uie complete model is illustrated in Figure A-5. 


A.3 NUCLEAR THERMOELECTRIC 

Nuclear primary power for electromagnetic propulsion requires 
thezmal-to-electrical converters capable of operating for long periods of 
time at high temperatures. Thermoelectric converters are in a unique posi¬ 
tion to meet these requirements. They are a direct energy conversion 
system. Static conversion and modularity are attractive characteristics. 
A large nember of parallel converters provide a system with high redundancy 
and reliability. Nodular design permits flexibility in design and 
operation. 

Mission studies have shown that a 400 kW nuclear thermoelectric power 
system can be deployed in a single shuttle launch with the largest payload 
of interest to this study. Thermoelectric power systems capable of 120 kW 
have been analyzed by JPL and LASL with the 400 kW being their ultimate 
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A specific mass model was generated from data in references 22 thru 32 . 
The referenced data was for different reactors, radiators, and operating 
strategies. The type of reactor and radiator has a strong impact on the 
mass of the system. The specific mass model was generated by curve-fitting 
the conservative data. The specific mass, for power levels ranging from 
1 kW to ISO kW (limitation of the available data) can be expressed as: 


«,p - 0.021 + 370.3/Ppp 


(A-3) 


where m^^ is in kg/V and Ppp in W. We will assxaie that the specific mass 
model (Equation A-3} is valid for power levels above the range of the 
available data. Equation A-3 and the referenced data points are presented 
in Figure A-6. 


^^Kelly, C. E., Ambrose, G. R,, 'Testing of the GPHS Electrically Heated 
Thermoelectric Converter,' lECEC, 82 9234, 1982. 

^^Eisner, N. B., 'Development of An Advanced RTG Using Segmented Thermoelec¬ 
trics,' lECEC, 1980. 

^^Koeing, D. R,, 'Heat Pipe Reactors for Space Power Applications,' AIAA-77- 
491, 1977. 

^^Ewell, R. C, Stapfer, G., 'Thermoelectric Conversion for Space Nuclear 
Power Systems,' lECEC, 82 9238, 19 82. 

^^Rockey, D. E., et al., 'Comparison of Evolving Photovoltaic and Nuclear 
Power Systems,' lECEC, 82 9011, 1982, 

^®Mahofkey, T., 'Future Space Power—The D.O.D. Perspective,' lECEC, 809016, 
' 980 . 

^^Ranken. W. A., 'Experimental Results for Space Nuclear Power Plant 
Design,' lECEC, 809142, 1980. 

^^Bnden, D., 'Space Nuclear Reactor Power Plants,' Los Alamos, LA8823-HS, 
19 80. 

^^Buden, D., 'Reactor Technology for Space Electrical Power (10 - 100 kW^)/ 
LA6891-MS. 1977. 
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Figure A-8. Specific Hass of Nuclear Brayton Cycle Power Systems. 

No qnaatitative specific mass data of the nuclear Rankine system was 
found in the literature. The mass of a nuclear Rankine system should be 
similar to the mass of the Brsyton system. Ihe Rankine cycle operates at 
higher temperature than Brayton cycle but there are no other significant 
differences. We assumed that the lower bound of the data for Brayton system 
is the mean for Rankine system. 

The specific mass is then expressed as: 

■,p - 0.013 + 178.6/Ppp (A-7) 

and is illustrated in Figure A-9. 

A.6 NUaSAR MHD 

An MHD energy converter consists of a gas dynamic nozzle (expander; 
followed by a duct. The duct is surrounded by a coil which produces a 
magnetic field in the duct. Electrodes are placed at the top and bottom of 






A.5 NUCLEAR BRAYTON/RANKINE 


TEe Brayton and Rankine cycles used in tnx bo-genera tor s have : ov spe¬ 
cific mass. ILey are designed as modular nnits and are connected in part, 
lei for redundancy and increased system reliability. Each nnit is asoal.y 
designed to deliver a fraction of the total pover required vith the capabil¬ 
ity of delivering much longer power in case of module failure. The input 
pressure level to either cycle is approximately proportional to the output 
power level required. 

Mission studies requiring power levels ranging from 25 to 600 kW have 
been conducted with the nuclear Brayton cycl e.^"^The power was 
delivered at 120 VDC or at 67 VAC. three phase. 1733 Hr. Reference 37 is a 
study of a Brayton cycle delivering 400 kW. SOO V, single phase, 3000 Hz. 
The voltage in all cases is limited by the alternator. As alternator 
technology improves, higher voltages and frequencies may be feasible. 

Data for the nuclear Brayton power system from references 32, 35, 37, 
38, 39, and 40 was used to model the specific mass. The specific mass model 


»sp = 0.019 + 190.6/Ppp 

for power levels above 1.5 kW and is presented in Figure A-8. 


(A-6) 


■^^Preston, J., 'React or / Bray ton Power Systems for Nuclear Electric Space¬ 
craft,' Princeton. 

3 8 

Layton, J. P., 'React or/Bray ton Power Systems for Nuclear Electric Space¬ 
craft,' lECEC, 809140, 1980. 

^^Gable, R, D., McCormick, J, E., 'Brayton Isotope Power System —The Versa¬ 
tile Dynamic Power Converter,' lECEC, 789383, 1978. 

^^Thonpson, R. E., 'Gas Cooled Reactors for Space Power Plants,' AIAA-77- 
490, 1977. 
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APPENDIX B 

FERFOSHANCE M0I£L1NG APKOACE 


B.l DnXODDCriON 

The thruster input pulse puraaeters sre (1) pulse energy, (2) cur¬ 

rent, J, (3) voltage, V, and (4) pnlae duration, tp. These input paraaeters 
are related to the thruster output paraaeters and aany other thruster char¬ 
acteristics. Ue thruster output paraaeters which affect the input require¬ 
ments are pulse aass and exhaust velocity. The other output paraaeters, 
pulse frequency and thruster lifetiae, do not directly affect thruster input 
pulse requirements. Other thruster characteristics (such as size, electrode 
geometry and auterial, propellant aaterial. etc.) clearly affect thruster 
performance and therefore must be considered. These other characteristics 
cannot be selected arbitrarily since they do affect thruster performance 
(i.e., thruster input requirements). In principle, they would be deter¬ 
mined by optimizing the design of the thruster. 

The efficiency of the thruster has a direct effect on the power and 
energy requirements, and therefore the size of the primary power supply and 
the ES/PC system. Therefore, the optiaization criterion for the thruster is 
the maximization of efficiency. Thruster characteristics aust be chosen to 
maxiaize efficiency. Assnaing this can be accoaplished, the input pulse 
paraaeters for optimua thrusters are functions of only the pulse aass and 
exhaust velocity. (That is, for each coabination of pulse aass and exhaust 
velocity there is only one thruster configuration which has aaximua effi¬ 
ciency and it is this one which is of interest here.) 

Let us assiae that a large body of data existed (either experimental or 
analytical) for optimized thrusters and that this data contained values for 
the thruster input pulse paraaeters over wide ranges of m^ and v^. A 
siaple, snilti-variate curve fitting technique could then be eaployed to 
develop the required relationships for E^, J, V, and t^ as functions of ap 
and v^. This is the approach which was used for modeling the electric rail 
gun.^ Table B-1 lists the thruster characteristics which were chosen opti¬ 
mally, leaving only m^ and v^ as independent thruster output paraaeters. 

^Bauer, D. P., Barber, J. P,, Vahlberg, C. J., 'The Electric Rail Gun for 
Space Propulsion,' NASA-Ql-165312, Final Report, February 1981. 
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TABLE B-1. 

PARAMETERS AFFECTING RAIL GDN EFFiaENCT 


rARviETER 


Increasing Effi¬ 
ciency Results 
'ROM Increasing/ 
Decreasing Mag¬ 
nitude of Para¬ 
meter 

Optimizing 
Value or 
Physically 
Imposed 
Constraint 

Physical 

Implications 

Value of 

Parameters 

USED OR 

Derived 


Rail Gun 

11 NDUCTANCE 

I. 'noth 


Increasing 

OPTIMim 

Rail guns with square 
bore are optimal.* 

L' • 0.63 mH/m 


length^ 

1 


Increasing 

Optimum 

Practical optimum 
length, X, is a func¬ 
tion OF BORE SI2E AND 
exhaust velocity for 
GIVEN values OF OTHER 
parameters. * 

X « kS.B X 


!rail Te.— 
iPESATuRE 


Decreasing 

Constrained 

■Tradeoffs required, 

TO OBTAIN PRACTICAL 

haste meat radiator 

SIZE AND ADEQUATE 

strength rails with 
acceptable res I ST IV- 

T - R50*K 



Specified by mission 
REQUIREMENTS. BOUNDS 
WERE IMPOSED. 


Accelerating stress 

LIMITED TO MAXIMUM 

allowable stress in 

PELLET.••• 

Both pellet density, 

F, AMD PELLET length. 
D, SHOULD BE MINI¬ 
MIZED. Practical 

pellet density CHOSEN. 

Pellet length in di¬ 
rection OF TRAVEL 
must BE LONGER THAN 
one-half the bore 

SIZE, TO PREVENT IN- 
BORE PELLET TUMBLING.* 


» SQUARE BORE AND PELLET LENGTH EQUAL TO ONE-HALF THE BORE S 
TO BORE SIZE (I.E., M • pH^/2 ■ 1500 h^). GIVEN THE PELLET 
INTERCHANGEABLY. 


lED optimal 


inductance. 











A tail gun sioulation conpnter program (containing these optimum 
characteristics) was nsed to generate a large data base of thruster input 
pulse requirements over wide ranges of m^ and v^. The data was then reduced 
using regression techniques to produce a set of parametric equations 
describing optimum EBG performance. The equations are: 

n « (B-1) 

J - 2.2xl0^mJ/3 (B-2) 

tp - 4.45xlO"®m^/^Ve (B-3) 

This simple data-based approach cannot be applied to the other three 
thruster types (PIT. TPP. MPO) because not enough perfonaance data (either 
experimental or analytical) is available over a wide enough range of m^ and 
Vg. In addition, it is not always clear that existing data represents 
optimum thruster performance. For example, data over considerable ranges of 
m^ and v^ i< usually generated experimentally with one thruster configura¬ 
tion. Intuitively, one woxild expect that the configuration of an optimn 
thruster would be different at these wide extremes. The task of generating 
a more complete body of performance data (either experimentally or analyt¬ 
ically) is considerable and is beyond the scope of this report. Methods of 
performance extrapolation based on a sound physical understanding of basic 
thruster operation most therefore be employed. 

In the following paragraphs we will examine each of the thruster input 
parameters and develop methods of extrapolating thruster performance based 
on existing performance data and physically reasonable assnptions concern¬ 
ing thruster operation. 

B.2 PULSE B4EBGT 

The purpose of this section is to relate the pulse energy to the 
independent output parameters of exhaust velocity and pulse mass. We begin 
with an assumption concerning the components of pulse energy. 
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Edis = /RJ^dt 


(B-7) 


where R is the apparent 'resistance* of the thrnster and J is the pulse 
current. 

For the BiPD thruster, the dissipative energy loss term was broken into 
two conponents; a resistive loss term and an 'electrode fall' loss term. 
The complete dissipative loss term can be written as: 

Edis ^ /EJ^dt + /VfJdt (B-8) 

where is the 'electrode fall voltage'. The first component is identical 
to the energy loss mechanism for the TPP and PIT. The second component is 
an energy loss attributed to 'electrode fall.' The voltage drop. is the 

voltage penalty that must be paid to remove current from the electrodes and 
is assumed to be constant. 

Thruster performance data strongly indicates that the predominant re¬ 
sistive element in the circuit is the plasma. We conclude that: 

Asaamptioa No. 3 - The resistive dissipation in the 
thruster is primarily due to plasma resistance. 

We further note that ionization must occur early in the acceleration 
(neutrals are accelerated only poorly at best). This implies that the final 
state of the plasma is reached in a small fraction of the pulse time and 
that: 


Aasvmptiom No. 4 - The apparent dissipative resis¬ 
tance remains constant during acceleration. 

Equation B-7 can then be rewritten as: 

Edis ' R/J^dt (B-9) 

(It should be noted that even if assumption no. 4 is not strictly obeyed, 
but the variation of resistance scales similarly with time as thrnster 
output is varied, then Equation B-9 will still be valid.) 
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Eqution B-8 can also be rewritten as: 


®dis * R/J^dt + Vf/jdt (B-10) 

Experimental thruster perfonance measurements indicate that for ex¬ 
haust velocities of interest 0 10,000 m/s), the thrust due to electrothermal 
effects is very small and: 

Assumption No. 5 - ^e acceleration of the plasma 
is primarily magnetic. 

Ihe electromagnetic impulse can then be written: 

■pV^ - (L’/2) / J^dt (B-11) 

where L' is the axial inductance gradient for the device. The inductance 
gradient is a constant which depends on device geometry but does not vary 
with time. 

Ihe dissipative energy loss term for the TPP and PIT can be rewritten 
by combining equations B-9 and B-11 thus producing: 

Edis “ (2R/L')mpV^ (B-l2) 

For the MPD, the resulting equation is: 

E^ij - (2R/L’) UpVg + Vj(2ipVp/L')^^^tp (B-13) 

These equations imply the very interesting result that the dissipated energy 
depends on the impulse bit. "pV^, 

We must now examine how R and L' vary with pulse mass and exhaust 
velocity. The plasma resistance should be independent of exhaust velocity 
according to assumptions 2 and 3 and will depend on plasma geometry and 
size. The resistance can be expressed as: 

E - ?£/A (B-14) 
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where ( is the plasas resistivity, Z is the effective plasms length, and A 
is the cross-sectional area of the plasma perpendicular to the path of 
current flow. 

The pnlse mass for the PIT and TPP is related to the plasma dimensions 
and plasma density by: 


Up - p.£A (B-15) 

wh^re t is the plasma length and A is the plasma cross-sectional area. 
Obvionsly, the plasma resistance, R, is related to the plasma mass, m^, 
through the geometry of the plasma and the material properties $ and p. 

We do not know, a priori, very much about the plasma geometry; however, 
it appears to be reasonable to assne that: 

Aaammptiom No. 6 - The plasma body will scale 

geometrically with pulse mass (i.e., ratios of 
dimensions roeain fixed as the thruster size and 
pnlse mass are varied). 

Assimiptions 2 and 4 imply that both plasma resistivity and plasma density 
are constant and independent of pulse mass and exhaust velocity. From 
Equation B-IS, we therefore conclude that plasma (and thruster) dimensions, 
d. scale as: 

d o m|/3 (B-16) 

and from Equation B-14 plasma resistance scales as: 

R o mp^^3 (B-17) 


which may be rewritten as: 


R . km 


-1/3 

P 


(B-lS) 


where k is a constant to be determined from experimental measurements. 
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Plasma dimensions and resistance scale somewhat differently for the HPD 
which operates in a qnasi-steady rather than a pulsed mode. 'Die MPD channel 
is completely filled with plasma. Die density and current have a stable 
distribution during the entire pulse. 'Die current flows radially in the 
plasma. From assnsptions 2, 3 and 4 we can write the plasma resistance in 

an MPD device in terms of an effective (constant) resistivity as: 


R 






o 

dr/(2jtr£) 

'i 


(D-19) 


where t ~ the length of the MPD channel 
r^ « inner electrode radius 
rg * outer electrode radius 

'Die integration may be completed to yield: 


a - ?la(rQ/r.)/2n.e (D-20) 

From assimption no. 6 we see that remains constant as the 

thruster is scaled and we conclude that: 

a o 1/ £ (B-21) 


From the data available on MPD thrusters we find that the mass flow rate 
scales as the annular area of the device: 

ip a rJ(l-(4./rp)^) (B-22) 

From assmption 6 we can therefore infer that: 

£ a r^ a r^ o (B-23) 

which may be combined with Equation (6-21) to yield: 

R - taap^/^ (B-24) 


where h is a com 


Lt to be determined from MPD performance data. 









Uie scaling relation for the inductance gradient, L', is the remaining 
nnhnovn. Ihe inductance gradient is a function of the geometry of the 
thruster electrodes. Parallel plate devices (such as the PIT) have a 
gradient given approximately by: 

L* - p^jw/h (B-25) 

where Pg ° permeability of free space 
w = separation of plates 
h « height of plates 

Ihe gradient for coaxial devices (such as the MFD) is: 


L' = (p^/2n)ln(r^/rj|^) (B-26) 

where rg ^ radius of outer electrode 
s radius of inner electrode 

From asstmption no. 6, L* is independent of pulse mass (thruster size) and 
exhaust velocity. 

Combining the inductance and resistance results for the PIT and TPP and 
collecting all of the constants together we obtain, (Equations B-12, B-18): 


^dis “ (B-27) 

where h is a constant which will be determined from performance data. For 
the MPD device we obtain a slightly different result, (Equations B-13, B- 
24) : 


Edis • 

where and k 2 constants to be determined from performance data. 

The pulse energy equation can now be rewritten for the PIT and TPP as: 


■ 


(B-29) 


and the efficiency as: 




















M 'CVS) 

Figaro B-3 . Efficiency Data and Model Resolts for tlie MPD. 

(1) Efficiency increases witb ezhaast Telocity (bat can never exceed 
lOO^l) . 

(2) The increase of efficiency with ezbanst velocity is most pro¬ 
nounced at low velocities (<10,000 m/s). 

(3) Efficiency increases with poise mass (thmstez size). 

(4) The effect of poise mass is pronounced below a critical size where 
efficiency drops rapidly with decreasing pulse mass. 

(5) At poise masses above the critical value, efficiency is insensi¬ 
tive to poise mass. 

B.3 CURRENT/POLSE TIME 

An average current can be defined from: 



Personal Common: 


Eenn Clark - Prince 


Dniversi ty. 







which can be further reduced by, substitution of Equation B-11, to: 

Javg = (2«pV,/L-tp)l/2 (B-34) 

Equation B-34 eaibodies most of the desired model characteristics. The 
inductance gradient. L', is independent of and y^. The average current 
therefore depends directly on the impulse bit. However, the average current 
and the pulse time are obviously interdependent. Assessment of current 
requirements must be preceded by an analysis of pulse time. Each of the 
thrusters (PIT. TPP and MPD) have different pulse time constraints and must 
be treated separately. 

B.3.1 The PIT 

The PIT accelerates a flat sheet of plasma by repulsion from a flat 
coil. As the separation between the plasma sheet and the coil increases, 
the mutual inductance (the PIT equivalent of L’) decreases. Ihe accelerat¬ 
ing force then drops and the plasma becomes 'decoupled* (i. e., current 
flowing in the drive coil no longer accelerates the plasma). Ihe data in 
the literature indicates relatively good (and constant) coupling up to a 
separation of about 1/10 of the coil diameter (the decoupling distance 
should scale directly with the coil diameter). Coil diameter will scale 
with pulse mass according to assumption no. 6 as expressed in Equation B-16. 
Correlating the available data^ we find that coil diameter, D, is: 

D » 66.5mi/3 (B-3 5) 

Therefore, the acceleration length, t, is given by: 

£ = 6.5ml/3 (B-36) 

The measured inductance gradient for the PIT is: 

L’ = 7nH/m (B-37) 


^Personal Communication - Dailey, C.L., ISW Space and Technology Group. 
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The accelera 


oltegi 


derived from Equation B-56 and is: 


- 2.07x10 “♦il/2^^3/2 (B-70) 

Hiese voltage ccaponents are shorn in Figure B-16. 



Exnoust velocity ikw/si 

Figure B-16. Voltage Components for the MPD. 


B.5 CGNCLOSIONS 

A number of specific conclusions can be drum from this investigstion 
of performance modeling for electromagnetic thrusters. 

(1) Thruster modeling for integrated propulsion system optimization 
requires that thruster input requirements (pulse energy, current, 
voltage, and time) be related to the primary propulsion system 
outputs (pulse mass, exhaust velocity, frequency, and total 
thrusting time). 

(2) Models which meet these requirements were successfully developed 
for all of the thrusters of interest. 

(3) The models developed for the thrusters (except for the EBG) were 
based on an understanding of the pl^sical principles of operation 
and a nimiber of asstmptions concerning the performance of opti¬ 
mized thrusters. 







B.4.3 The MFD 






of tlie total poise tiae. Hie ioaizatioo voltage can then be written froa 
Eqoa tion B-53 as: 

Vion ' 7.5 x10®»J/3/v^ (B-64) 

Ihe acceleration voltage is derived froa Eqoation B-5S as: 

Vacc = 0.311aJ/3v^ (B-65) 

The dissipation voltage, coapoted frca Equation B-57 is a constant: 

'^dis * *7.2 (B-66) 

These voltage coaponents are shown in Figure B-14. The TPP thruster does 
not require a 'breakdown voltage', since the breakdown across the electrodes 
is initiated by an external ignition circuit. There is a ainiana allowable 
voltage for operation of the device. Lack of experiaental data precluded 
developing a relationship for this paraaeter. 













Exnoust velocity 'ra/si 

Fignre B-11. Voltage Coaponentt for the PIT. 

' ionization voltage' but reflects the practical reqnireaent of achieving 
adequate electric field to ionize the gat. Fron the ezperiaental data. 
Equation B-53 was reevaluated to yield: 

Vbdn - 7.5*10^»>J^3 

which is plotted in Figure ft-12. Cosparison of Equation B-62 with Equation 
B-S9 shows that the breakdown voltage is higher than our 'ionization volt¬ 
age'. This indicates that ionization will be coapleted in lest than 10% of 
the pulse tiaie. 

B.4.2 The TPP 

The average voltage for the TPP thruster can be found froai Equation 
B^52 and is: 

V^g - 0.311m^^^^e * 7 .SxIo'^bJ/ 3/v^ + 87.2 (B-63) 


which is shown plotted in Figure B-13. 






Figoze B-10. Average Voltage for the PIT. 


Eacli of the voltage conpooeats was ezaiined separately. Assoaing that 
e plasata is foraed during the first 10% of the total poise duration, the 
ctor (tp/t) is equal to 10. Hie ionization voltage can be written fron 
na tion B-53 as: 


Vion « 4.61 x10^-J/3 (B-59) 

e acceleration voltage derived frea Equation B-55 is: 

V.gc - 2.59xl0~'*al/3^ (B-60) 

d the dissipation voltage derived fm Equation B-S7 is: 

Vdis - 7.52xlO”^Vg (B-61) 

e three voltage ccaponents are plotted in Figure B-11. 

The PIT also hat a ainiata allowable voltage required for gas breakdowj 
d plasaa fozaation. This 'breakdown voltage' is identical to oni 










Tbe ionizmtion voltage. as defined here. ahoold not be confused with the 


voltage required to ionize (the breakdown voltage). Ihe factor. (t/tp>i it 
the fraction of the pulse duration in which the ionization takes place. For 
example, the PIT and TPP require a very high voltage for plasma formation at 
the beginning of the pulse. Once the plasma has formed, this voltage is no 
longer required. On the other hand, the MPD thruster requires voltage for 
ionization of propellant throughout the pulse. For the NPD. (t^/t) > 1 and 
the 'ionization' voltage is written as: 

^ion * Vi^^ 

The voltage required to deliver the acceleration power to the pulse 
mass is the acceleration (or useful) voltage. It is related to the kinetic 
energy, the average current, and pulse duration as: 

'^acc = 

For the MPD the acceleration is quasi-steady and the acceleration voltage 
can be written as: 


Vacc « Ve/J 

The voltage that is associated with the energy dissipated during the 
pulse is related to the dissipated energy, the average current, and pulse 
duration. The dissipative resistance of the thrusters was assumed to be a 
constant throughout the pulse. The dissipated voltage is given by: 

V. - 

We now exsmine these voltages for each thruster. 

B.4.1 The PIT 

Since all of the factors in Equation B-S2 are known, the PIT average 
voltage can be written as; 

V — . - 2 . 59 xl 0 ”^mi/^v^ + 4 .dlxlO^m?;^^ -t- 7.52xl0~^v. (B-58) 







20 30 40 


Figure B-9. Current for tie MPD. 


The avenge voltage is conputed fron the pulse energy. the pulse time 
the avenge current as: 


origin of the voltage can be better understood if it is written as the 
of voltage components related to the components of pulse energy. These 
tage c(Mponents are denoted by: ^ion ~ ▼oltage required to heat and 

Lze the propellant, - voltage required to accelerate the propellant, 

~ the voltage characteristic of the energy dissipated. A minimna 
tage required to initiate operation (breakdown vol tage) has also beei 
lyzed for some devices. 

The 'ionization' voltage is computed from the ionization energy 
□ation (B-6)} substituted into Equation (B-52} to yield: 


n ei/(Jt„) z (t /t) 


(B-53) 








B.3.3 The MPD 


Since the MPD operates in a 'qnasi-steady' node, the pulse duration is 
an independent parameter. The current is therefore independent of the pulse 
time and is assmed constant during the pulse. The electromagnetic thrust 
for an MPD thruster can be written aa: 

F - « L’J^/2 (B-48) 

This can be generalized to form a relationship between the current and 
the independent pulse parameters: 

J a (B-49) 

or written in equation form as: 

where h is an experimentally determined constant. From the awailable data^ 
the constant was evaluated and Equation B-SO becomes: 

J - 2.42xl03(*j,Vg)^^^ (B-51) 

which is illustrated in Figure B-9. 


B.4 VOLTABE 

The thrusters treated are all electromagnetic and are 'current driven' 
(i.e., the acceleration force is related to current, not voltage). These 
devices do not, in general, require a high voltage to accelerate the pulse 
mass (although some require a relatively high voltage for the formation of 
the current-carrying plasma and to produce the required short pulses). The 
voltage is, therefore, a derived parameter in electromagnetic thrusters and 
is indicative of the impedance of the device. The primary function of the 
ES/PC system is to provide the required current. 


























Exnaust velocity (wn/s) 

Figare B-5. Average Corrent for the PIT. 


measured ptilse times are, as expected, longer than Equation B-'39 predicts, 
but only by about 15%. There are a few measurosents of peak current 
reported and, as expected, they correlate with pulse mass and exhaust 
velocity as predicted by Eqnation B-41. The data'^'S for peak current 
yiel ds; 

J„ax = 

which is illustrated in Figure B^. 

Comparison of Equations B-40 and B-41 indicate that, for the particular 
ES/PC system used in the experiments, the peak-to-aver age current ratio was 
approximately 1.75. 


^Dailey, C.L., 'Large Diameter Inductive Plasma Thrusters,’ AIAA-79-2093, 
Presented at the 14th International Electric Propulsion Conference (Prince¬ 
ton, NT) . 

^Dailey, C. L., 'Large Diameter Inductive Plasma Thruster,' AIAA-79-2093, 
presented at the 14th International Electric Propulsion Conference (Prince- 
t on, MI) . 
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SECTION C.l 


INTEOOOCTION 

Indactors can b« naad eitJiar as energy stores or as pal se-f orming 
coaponents in poise power systeas. In tbe energy storage role indactors are 
attractiwe becaase they offer high stored energy at relatively low aass. In 
addition, they provide a 1 ow-iapedance, constant-oorrent sooxce characteris¬ 
tic. This characteristic is very iaportant for soae applications sach as 
el ectrcaagnetic laonchers and electroaagnetic thmsters. The disadvantages 
of indactors as energy stores is that they are lossy (i.e., have a low Q) 
and are therefore relatively inefficient. The losses can be eliainated and 
the efficiency iaproved by eaploying snpercondnetors. Snpercondactors, 
however, are tot tolerant to poise discharge and pose serioos thezaal design 
probleas in high corrent applications. 

Indactors sre also ased as pal se-foraing coaponents to control the 
shape and dnration of poises delivered froa capacitive energy stores. In¬ 
dactors can be osed in either a noraal LC coabination or in anlti-staged 
poise foxaing networks. In these applications, the indactors carry current 
for only a short tiae and the resistive losses inherent in indactors are not 
so important. 

In aany applications of interest for both inductive energy storage and 
inductive pulse foxaing, ccaponent aass most be ainiaized. The objective of 
this report is to develop methods of ainiav aass selection and sizing of 
indactors for both energy storage and pulse foxaing duty. 
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SECTION C.2 


CIRCOIT CHAEACmiSTICS 

Tha circuit perforausco of inductors pisys an important role 
tor sizing and selection. The energy efficiency of the induct 
circuit IS of particular interest. Ve define the efficiency as: 

n - ^use/^^dis^^use^ 

where E^^^, * dissipated energy. 

E^f, * useful energy extracted. 

13ie dissipated energy is: 

9 

Edis - B/I^dt (C-2) 

where R « resistance of the inductor. 

The useful energy is: I 

E^se - UlJ - l|)/2 (C-3) 

where Iq * mazimom current during period of interest p 

If = ainimna current during period of interest 
Substituting Equations 2 and 3 into Equation 1 prowides; 

n - 1/(1 + 20/x) (C-4) I 

where 0 » /l^dt/CIg - If), the efficiency factor 
X ^ UR. the inductor time constant. 

The factor, 0, depends on the details of the application and the character- 0 

istics of the circuit. 

In pul se-f oraing duty, 0 can be readily evaluated from the circuit 









In anergy storaga duty, 0 is related to the frequency of pulsing and 
the depth of discharge. Ihe basic current-tiae characteristic for energy 
storage duty is illustrated in Figure C-1. The discharge time is typically 
orders of magnitude shorter than the charging time. The efficiency of the 
inductor is therefore primarily determined by the charging process. TVo 
different charging yariations are illustrated, a linear current rise and a 
square root current rise (which more nearly approximates the desired con¬ 
stant power charging). The factor 9 can be eyaluated for both cases and is 
related to the frequency by: 


where a = an efficiency factor. The efficiency factor is a function of the 
current ratio, and is illustrated in Figure C-2. For reasonable 

current ratios (<0.5) the efficiency factor is independent of current ratio 
and is approximately equal to one. 
















SECTION C.3 


INDOCTOR SIZING 


C.3.1 THE INDDCTANa: 

Froi tEe point of viaw of indactoc sizing and design, the inductance is 
not a design paraaeter. TEe designer is not free to choose the inductance 
in order to niniaize aaaa. The inductance is usually fixed by the appli¬ 
cation and the designer's job is to size and design an inductor with the 
required inductance and ainiaiai aass. Frca the designer's point of view, 
therefore, the inductance is an independent paraaeter. 

The inductance of an inductor is dependent upon the configuration and 
the size of the inductor. The inductance may be conveniently expressed as: 

L » kNz (C-6) 


where L • inductance 

k » inductance factor 
N a nmber of turns 
z - length of conductor. 

The geoaetrical factor, k, depends on the choice of configuration (toroid, 
solenoid, etc.) and on the geometry of the inductor (length-to-radius ratio, 
etc.). The value of the inductance factor is independent of coil size. 
Values for the inductance factor can readily be calculated for all ccmaon 
configurations and geometries. In this report we concentrate on soleniods 
as they are simple and display the same basic scaling characteristics as 
other multi-turn concepts such as toroids. 

The geometry of a solenoid is illustrated in Figure C-4. The induc¬ 
tance tactor for solenoids can readily be evaluated from the formulas of 
Grover^ and Wheeler.^ The inductance factor is a lunction of the geometri¬ 
cal ratios c/a and b/a and is illustrated in Figure C-S. We note that for 


^Grover, F. W., 'Inductance Calculations,' 0. Van Nostrand Co., New York, 
1946. 

^Wheeler, H. A., 'Inductance Formulas for Circular and Square Coils,' IEEE 
Proc., 70, pp 144 9-145 0, 1982. 
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Figure C-4. Solenoid Geoaetry. 


long coils (b/s > 1) tbe inductsiica factor is not a strong function of tbe 
tbickness ratio (c/a). Ibe parameters for a Brooks coil are also identified j 

in Figure C-5. The Brooks coil is important as it is the solenoid config¬ 
uration which has the highest inductance for a given length of conductor. 

It is soeietimes thought of as a minimum mass configuration. However, it is 
only a minimum length configuration. As we shall see in the following 
sections, inductor mass depends on other factors as well as conductor 
leng th. 

The inductance factor for solenoids can be expressed approximately as: 

k > k^ - ki ( 

where k^ = 2xl0"‘^ [2.78/( (b/a)+l)+ln(l-K) .39(b/a)) ] 
kj^ - 6.2 8xl0“^(*/b)(c/a)k' 

k- = (0.331-0.011/(b/a))exp[-(0.361+0.088/(b/a))(c/a)]. 

The approximation is very good over the range of geometries of interest. 






Radi us/length Ratio (a/b) 


Figue C-S. ladncta&ce Factor for Solenoids. 

C.3.2 INOOCTOR MASS 

The Bass of a solenoid can be expressed as: 

a a pAz (C-8) 

where a > inductor aass 

p a onss density of the conductor 
A a cross-sectional area. 

Inductor aass is ainiaized by ainiaizing the product of z and A. Inspection 
of Equation C-6 shows that z will be ainiaized if N is aaxiaized (for a 
fixed inductance) and that inductor aass will be directly related to the 
inductance. Equations C-6 and C-8 prowide no information on the selection 
of conductor cross section, A. For a ainiaua aass inductor, A must be as 
smal 1 as possibl e. 

C.3.3 CONDDCrOS OIOSS-SECTIONAL AREA 

The conductor cross-sectional ares can in principle be reduced until 
the stress and/or the teaperatnre (due to resistive dissipation) becomes 
unacceptable. Further reduction in cross section will result in coil 
failure. 

C.3.3.1 Thermal Constraints 

The thermal design of an inductor can proceed in either of two direc¬ 
tions. In the first approach the inductor is not actively cooled, but is 
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small c/a (<0^}. Hie effect of thickness-to-radias ratio is illustrated 
more clearly in Figure C-9. 

C.4.3 CONDOCTOR TEMPERAIURE 

The conduotor temperature primarily affects the resistivity but also 
bas an effect on beat transfer. Hie variation of specific energy vi tb 
temperature is illustrated in Figure C-10. Hie specific energy of tbe 
inductor is almost independent of temperature. Hie total mass of tbe system 
must include tbe cryogenic system and will not necessarily show tbe same 
effects as tbe indnctor only. 



Temoerature iK) 


Figure C-10. Specific Energy vs. Conductor Temperat'jre for Al-1100-0. 


C.4.4 CURRENT 

Tbe effect of current on inductor specific energy is illustrated in 
Figure C-11. At low current, tbe specific energy is proportional to tbe 
current and bas a siigbt dependence on tbe inductance. At bigb current, tbe 
energy density is independent of current an 


id inductani 






Aspect Ratio 8/A 

Sp«ci£ic Energy vs. Length-to-Rndi 




C.4.1.3 Tlio Geometry 


Tile geometry was varied to determine inductor mass sensitivity to these 
factors. The ranges examined were 

0.05 < c/a < 1 (C-32) 

0.125 < b/a < 16 (C-33) 

s/A » 2500 m-1 (C-34) 

C.4.1.4 The Outputs 

The program calculates a number of inductor parameters. For the pur¬ 
poses of comparison the specific energy (stored energy per unit mass) was 
chosen as a figure of merit. The higher the specific energy the 'better' 
the inductor (i.e., a lower mass for a given stored energy). This parameter 
also permits rapid comparison of inductors with other types of energy stor¬ 
age for which specific energy is usually known. The specific energy was 
examined over the entire range of parameters described in the preceding 
paragraphs. 

C.4.2 INIXKH'Oa GEOMETRY 

The 1 ength-to-radius and the thicJmess-to-radius ratio both have an 
effect on the specific energy of inductors. A typical electromagnetic 
thruster case was chosen to examine the effect. The inductance was lOufl, 
the current was lOOhA, the conduction time was set at lO^s (an energy 
storage case) and a temperature of 300^K was assumed. The results are shown 
in Figures C-8 and C-9. 

Figure C-8 shows the effect of length to radius ratio at various thick¬ 
ness to radius ratios. For these conditions the specific energy is a few 
hundred joules per kilogram, an order of magnitude higher than dielectric 
capacitors. The specific energy for very long solenoids is approximately 
1/2 that of very short ones. Increasing the th i c kne s s-t o-r adi us ratio 
improves the specific energy, however the results for very short (b/a < 1) 
solenoids are not accurate. The formula currently used to compute the 
inductance factor, k, is only accurate for either large b/a (>1) and/or 






SECTION C.4 


RESULTS 

Hie compater progm wms exercised over s vide range of inductor param¬ 
eters. The purpose was to examine the ability of the program to select 
realistic minimum mass inductors and to determine the characteristics of 
such inductors. 


C.4.1 THE PARAMETERS 


C.4.1.1 The Independent Parameters 

The independent parameters are the inductance, L. the current, I, and 
the conduction time, t. The ranges investigated were derived from electro¬ 
magnetic thruster considerations and were: 


O.lMfl < L < 25|iH (C-29) 

loH < I < lO^A (C-3 0) 


10-«s < t < lOyr {C-31) 

C.4.1.2 The Material Properties 

The properties of Al-1100-0 aluminum were used. The performance at 
three temperatures, 21°X, 77°K and 300°K examined. The properties used 
are summarized in TABLE 1. 


TABLE 1. MATERIAL PROPERTIES. AL-1100-0 

Temperature Density Stress Resistivity Temperature Action Heat 

Rise Constant Transfer 

(OR) (kg/m2) (MN/m3) (0-m) (OR) (A2-s/m'‘) (w/m^) 


21.00 2700.00 70,0 O.llOE-8 50.0 0.130E17 0.100E5 

77.00 2700,00 70.0 0,400E-8 50.0 0.110E17 0.500E5 

300.00 2700.00 70.0 0.255E-7 50.0 0.300E16 0.100E7 
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C.3 .5 COMPUTATIONAL TECHNIQUE 


TEa aiaiatai aass configuration is found by sianltaneons solntion of the 
equations deTelopad in tEe preyicns sections. A closed fom solution is not 
possible and a nnaarical approach anst be used. A computer program was 
developed to find the solutions. lEe inductance, current, and conduction 
time vera treated as the independent variables. The material properties of 
resistivity, density, action constant, and maximum stress must be chosen. 
The thermal transfer rate and periaeter-to-area ratio must also be selected. 
Finally, the geometrical factors c/a and b/a must be chosen. Ihe program 
then finds the miniaia aass inductor and determines the mass, resistance, 
number of turns, stress, and dimensions of the inductor. The choice between 
active and passive cooling is made based on minimizing mass. Reduced resis¬ 
tance configurations are found by incrementing the conductor cross sectional 
area and repeating the computation. The stress and thermal requirements are 
reduced in reduced-resistance inductors. 
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i&dactance during rapid discharge. During charging, current and magnetic 
field diffuse into the conductors. When the inductor is subsei^uently dis¬ 
charged, the magnetic field diffused into the conductors (and the energy 
associated with it) cannot be recovered. Only the field and energy stored 
in the 'open* area of the inductor and in the skin depth of the conductors 
can be recovered. The result is an apparent drop of inductance during 
discharge. The effect is most pronouaced in thick coils, such as the Brooks 
configuration, where reductions of 20% or more can be expected. In addi¬ 
tion. for multi-layer coils, internal circulating currents are induced which 
result in significantly increased resistive losses. 

These high frequency effects can be controlled by minimizing inductor 
thickness, c, and by maximizing inductor thickness-to-radius ratio, c/a. If 
c/a is constrained to be less than about OJ., the high frequency effects are 
modest. 

C3.4 INDOCrOR RESISTANCE 
C^AJ. Minimum Mass Inductors 

In the preceding paragraphs the conductor geoaetry for minimum mass 
inductors was developed. The resistance of these inductors can be readily 
cal cul ated from 

i^^z/A (C-28) 

The cross sectional area. A, is deteained from thermal constraints and the 
conductor length, z, is determined from stress considerations. The minimum 
mass constraints result in maxima resistance. For some applications this 
resistance may be unacceptably high. The resistance can always be reduced 
but there is a mass penalty. 

C.1.4.2 Resistance-Mass Scaling 

In order to reduce the resistance the size of the inductor m’lst be 
increased. The process begins by increasing the cross sectional area. The 
size and other characteristics of the inductor can then be recalculated 
using the previously developed formulas. There is no simple relationship 
between inductor mass and resistance. 
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AI (1100-0) 



Figura C-7. Tha Electrical Skia Depth for Alomizitia. 

CJJ.4.1 Resistaaoa 

For poises ia which the skin depth is less than the coadoctor thick¬ 
ness, the resistance of the coadoctor increases sbowe the DC resistance 
inyers^.'r with the skin depth. Iherefors, an indnetor can have a low DC 
resistance hot 7et have a high discharge resistance. 

In energy storage ' appl ica tions the indoctor most respond in two quite 
different time scales. Charging is a relatively slow process related to the 
frequency at which poises are reqoired. The discharge time, or poise dura¬ 
tion, is uaoally ooch shorter (typically by orders of aagnitode) than the 
charging tiae. (Otherwise an energy store would not be required.) The skin 
effect resistance is not usually a problea during charging, but can be 
iaportant during discharge. The thinner the conductors are, the less depen¬ 
dent the resistance is on pulse tiae. The high frequency resistance of an 
indnetor can usnal ly be designed to acceptably low values without seriously 
affecting indnetor sizing or configuration. 

C.3.3.4.2 Inductance 

The other effect of diffusion 1 iaited skin depth is a reduction of 
inductance. This effect can result in a substantial reduction of -apparent 
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forces the coil thickness, c, to increase (in order to retain the required 
conductor cross-sectional area). Under these circumstances, the required 
inductor thickness, c, can be considerably greater than the mean inductor 
radius, a. This of course is an impossible configuration. The stress (and 
therefore, the number of turns) must be reduced in order to achieve a 
geometrically possible inductor. The implication is that for certain appli¬ 
cation regimes a highly stressed inductor cannot physically be constructed. 
The geometrical constraint may be expressed as: 

c/a < 1 (C-26) 


and must always be applied. 

C.3.3.4 High Frequency Effects 

Both energy-St or age and pulse-forming inductors are subjected to high 
frequencies during pulsing. The effects of high frequencies on the resis¬ 
tance and the inductance of the coil can be important, particularly for high 
current designs where conductor size is large. When current is changed 
rapidly in a conductor, current diffusion effects confine the changing 
current to a thin surface layer or skin. The depth of this layer is related 
to the pulse duration by: 

Cj = (C-27) 


where c, =* »hin depth 
fp * pulse time. 

The skin depth for aluminum at various temperatures is illustrated in Figure 
C-7. The skin depth for millisecond pulses (those of interest here) is 
less than 1 cat. If the skin depth is less than the conductor thickness then 
both the resistance and the inductance of the coil are affected. 


C-13 




















The presence of fins and/or cooling channels increases the geometrical 
cross section of the conductors. The geometrical cross section and the 
physical cross section are then related by: 

A =■ pAg (C-18) 

where Ag ° geometrical cross section 
P = 'solidity' ratio ^1.0. 

The solidity ratio could be as low as 0.5 in order to achieve mazifflum 
cooling rates. 

C.3.3.2 Stress Constraints 

A solenoid acts as a magnetic pressure vessel. The magnetic field 
exerts a bursting pressure on the solenoid and this in tnrn induces stress 
in the windings. The stress in the conductors of a solenoid is related to 
the magnetic pressure by: 


a - Paa/cp 


(C-19) 


where « = conductor stress 

Pg = magnetic pressure. 

This equation is strictly true only if the stress is uniform throughout the 
conductors. The magnetic pressure is related to the ma6>^etic flux density 
by: 


Pa = (C-20) 

where B = magnetic flux density 

= permeability of free space. 

The magnetic flux density is related to inductor geometry and current 
through the Biot-Savart law by: 


B = UoNI/b (C-21) 

Combining Equations C-19, C-20, and C-21 we obtain: 
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by increasing tbe beat transfer rate (i.e., forced convection/nncl ear boil¬ 
ing), and by increasing tbe perimater-to-area ratio (i.e., finning or cbanr- 
neling of tbe conductors). 

Heat transfer rates Tary froa a few kf/n? for cryogenic cooling to a 
few MV/b^ for boiling water cooling. Xbe aaxian acbiewable transfer rates 
depend strongly on tbe detailed design of tbe coolant channels and tbe 
▼ elocity, pressure and Reynolds Nneber of tbs coolant flow. We shall assnae 
tbe following beat transfer walnes for this study: 

4h - loS/a? (C-13) 

for liquid lydrogen, 

4jl - 5xl0'*w/a2 (C-14) 


for liquid nitrogen, and 


q^ - lO^w/a? (C-15) 

for water. 

Ibe periaeter-to-area ratio, s/A, is maximized by finning and/or by 
providing many small channels for coolant flow in tbe conductor. In addi¬ 
tion tbe surface of tbe fins or channels should be finely grooved to farther 
increase tbe available cooling area. If this is done, tbe ratio is related 
only to tbe mean thickness of tbe conductor fins or webs by: 

s/A - 4/d (C-16) 

where d > mean thickness of tbe fins or webs. Xbe perimeter to area ratio 
is maximized by minimizing tbe fin or web thickness. Xbe minimom thickness 
is in practice probably related to the mean conductor thickness (i.e., c in 
Figure C-4). We shall assume that tbe minimum is given by: 

d - 0.1c (C-17) 
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Figure C-6. Action Constant for Alminim. 

C.3.3.1.2 Active Cooling 

Active cooling of a conductor can be characterized by equating 
resistive dissipated power to the thezaal transfer fron the cooled surfa 
This nay be written as: 

1^5z/A - q«z (C- 

vhere { - resistivity 

q > heat transfer rate 

s > 'wetted' perimeter available for cooling. 

Equation C-11 can be rewritten as: 

A = I(5/4(s/A))^''2 (C- 

Ihe perineter-to-area ratio, s/A, is a function of the geometry of 
conductors. The cr os s-se ct i ona 1 area (and therefore, inductor mass) 


the 
ce s. 

- 11 ) 


-12) 

the 
^ is 
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designed as a themal inertia. Ilie conductor temperature is allowed to rise 
and the cross section is selected to keep the temperature rise acceptable 
for the required current level and the required conduction time. In the 
second approach, the conductors are actively cooled and the cross section is 
sized to permit steady-state transfer of the dissipated energy into the 
coolant. Ihe approach which vil 1 provide the minimum-mass coil depends very 
strongly on the application. The tradeoffs and proper selection are made 
during computation. 

C.3JJ.J. Thermal Inertia 

The required cross-sectional area for a thermal inertia design can be 
readily determined from the action integral: 

/I^dt - gjA^ (C-9) 


where I > current 

t > time of conduction 

• action constant (dependent on temperature rise) 

A a cross-sectional area. 

The action constant, g^^, is a function of the conductor material, the ini¬ 
tial temperature, and the desired temperature rise. The values for aluminum 
are illustrated in Figure C-€. 

Assuming a constant current. Equation C-9 can be rewritten as: 

A>I(t/gi)l^^ (C-10) 

The cross-sectional area, and therefore, the inductor mass, can be reduced 
by reducing the current or the pulse time or by increasing the action 
constant. The action constant can be increased by material selection, by 
increasing the allowable temperature rise, or by decreasing the initial 
temperature. Copper has a higher action constant than aluminunr, however, 
that is more than compensated by the higher density of copper. Aluminum is 
usually the best overall material choice for minirntmi mass inductors. 
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C.4.6 CCNDOCriQN TIME 


ni« specific 9ti0rgy is also a foactioa of the conduction tine and is 
shown in Figure C'13. In pal se-f oning duty, tiae corresponds to the poise 
tine. In energy storage application, condnction is nsoally continuoas or 
nearly so. For very short pulse tiaea the inductor is sized as a theraal 
inertia. The ainiaua conductor cross section is then directly related to 
the pulse tiae. Ihe specific energy of such inductors can be yery high. As 


pulse tiae is increased the inductor auat be actiyely cooled and the conduc¬ 
tor cross section is independent of condnction tiae. The specific energy is 
usually lower for long-pulse actiyely-cool ed inductors than for short-pulse 
thenaal inertia designs. 










SECTION C.5 


OONaOSIONS 


S«v«zal impoztaat coaclnaiona can be drawn fm this study of indnc*- 

tors. The conclusions aze: 

(1) The specific energy of noraal inductors is very attrsctiwe when coa- 
pazed to that of dielectric capacitors and is coaparable to inertial 
stores. 

(2) Cryogenic operation has little effect on the specific energy of induc¬ 
tive stores. Uie aass of the cryogenic systea anst also be included to 
arrive at a systea level conclusion. Cryogenic operation does, how¬ 
ever, significantly rednce the inductor resistance. 

(3) The aass of an inductive store is a very coaplex function of material, 
operational, and application paraaeters. It is not possible to arrive 
at a single paraaeter (such as specific energy) which will completely 
characterixe inductive energy stores for any application. The aass 
must be evaluated specifically for each application. The coapntational 
technique developed in this study meets this need. 
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APPENDIX D 


PARAMETRIC MODEL OF SUFERCONODCTING ENERGY STORE 
FOR SPACE PROPULSION 


Sapercondacting inductiTO eiMXgy storage is an alternative method for 
long-term energy storage with the inherent advantage that the energy is 
immediately available in electrical form from a low impedance source. Com¬ 
pared to storage inductors made from normally conducting metals, supercon¬ 
ducting inductors have no joule heating and an essentially infinite time 
constant. However, there are losses that vary with the rate of field change 
in the superconductor and there are refrigerator power demands that depend 


sn these losses, on the surface area of the dewar, and on the currents that 
Bust flow fron ambient to cryogenic temperatures. Maximum current densities 
and optimal field strengths are determined by the critical parameters of the 
superconductor. For large high-field systems, the mass and size of the 
^orce containment structure become the limiting factors. Temperature may be 
dictated by operational considerations or by refrigeration costs. 

The objective of this note is to develop a simple parametric model of a 
superconducting energy store for space propulsion systems so as to permit a 


txons of a Kilojoule to half a megajoule, currents between IS and 1200 KA 
roltages generally below 1 kV, pulse lengths of 1 (is to 6 ms, and pul si 
repetition frequencies between 0.4 and 19 Hz. 

This comparison is somewhat complicated by the fact that two differen 
types of energy efficiencies can be defined for a cryogenic system: 

(1) An overall efficiency, based on total power consumption 

including that of the refrigerator, and 

(2) A conversion efficiency. ^2* which only looks at the 'higl 
quality' pulse power inputs and outputs and disregards the powe 
consumption of the refrigerator or the energy content of any 
'used' ciryogenic fluid. 


The efficiency ^2 will be more or le 


for a given oper 


the repetition rate decreases and when long period 








D.1.1 Coarooneats of SnuTcondacting Baargy Store 


Tii9 aaia coapoaaats of a saparcoadactiag saaxgj storage system are: 

1. Hie sapexcoadaetiag ladaetor, 

2. Hie earreat leads iato tlie oxyogeaie regioa, 

3. The dewar systea, aad 

4. Hie refrigerator. 


D.2 THB SOmCamOCTINB HOIDCIOE 
0.2.1 Oaeratiat Coaaideratioaa 

To fully specify a saparcoadactiag storage systea that vill deliver a 
specified eaergy palse at a specified repetitioa rate, requires a aiaber of 
iteratioaa oa the folloviag variaolea or paraaeters: 

OptiauB field straagth for aiaiaia size aad suss, 

° Field shape aad friagiag as dictated by operatioaal coasideratioas, 

^ Ratio of stared to delivered eaergy, 

^ Type of sx^ercoadactor aad its operatiag teaperatare, 

0 Fora aad saoaat of stabilizer, 

** lateraal desiga of the ooadactor to aiaiaize hysteresis aad eddy 
earreat losses, aad 

° Coolaat flow, plus locatioa aad siziag of the cool iag chaaaels. 

Clearly, this is soaevhat beyoad the scope of tae preseat analysis. 
Hovever, it is always true that the highest possible field will give the 
aost coBpact storage. Nevertheless, wheaever high carreats aad high repeti¬ 
tioa rates are required ia a sel f-saf ficieat env ironeat, ref rigeratioa 
ecoacaics, as deteraiaad by the Carnot factor, will pash the systea to the 
highest possible teaperatare. 

Most sapercoadactiag aagnets built to date operate at a teaperatare of 
about 4.2 C ia liquid heliaa boiliag at ataospheric pressure. Superconduc- 
tiag systeas with dedicated closed cycle refrigerators perait operatioa at 
any teaperatare between 2 and 20 K. For the present task, an operatiag 
teaperatare of 13 C is suggested, ia conjaaction with a refrigerator ontpnt 
at 12 I, This is not a conventional choice for either aagnets or storage 
indactors, bat it is feasible as Stevenson (1973)^ has reported a S T gas 
cooled solenoid operating at 13 E. 
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D.2.2 M««« of Inductor 


In tho following Oiscnssion tlio nsunl notation for solenoidal coils 
will bo used. Hie ratio of outer to inner radins is 

“ “ *2 / • ( 2 - 1 ) 
and the ratio of half-length to inner radins 

P * b / Sj, . (2-2) 

Very large superconducting energy storage systems have been proposed 
over the past IS years by Roger Boos and his associates at the University of 
Wisconsin. The latest of these reported by Eyssa, Boon, et al.2 is a 
100 hWh (360 Ml) study for NASA. It uses an earlier parametric analysis by 
Moses^ and suggests 10 T. 2 C coils at current densities of 2S0 to 
500 A/nm2^ also accepts the difficult task of noit-uniform distribution 

of superconductor, as dictated by local field strengths, and proposes novel 
puce tension structures for force support. Since the design is not con¬ 
cerned about stray magnetic fields, it favors very short solenoids 
(P « OJ), where most of the stored energy is outside the winding. 

For smaller systems, one can afford tj be somewhat less optimistic and 
consider longer coils with uniform current distributions. The mass of the 
inductor can be divided into three conponents: the superconductor, the 

stabilizer, and the structure providing force containment, support, and flow 
ch’!*<nels for the cryogen. 

Realistic asstaptions would be that the superconductor proper, with a 
density of 8 gm/cm^ will pexmit a current density, J^, of 10^ A/ii|2 (or 1000 
A/mn^), and that the stabilizer and channel structure of density, d^, will 
be present at an 'equivalent current density'. J, of 1 to S'lO^ A/ib 2 (lOO to 
500 A/sa?). The value of d^ would be 9000 kg/m^ for copper and stainless 
steel, and 2700 kg/m^ for alminum and its alloys. 

A still simpler approach can be based on an analysis of Label!,'( who 
shows that for solenoids with uniform current density the specific mass. 
M/Ej, is proportional to the average density, d^^^, (usually close to d^) 
divided by the square of the maximum field on the conductor, times some 
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(grsplied) coefficients tiist depend on a end 0. For sasll magnets (Ej< qj 
MJ), the mass of the winding and coil mandrel predominates, and at d^ * 9000 
kg/m3 


/ (E,*B^ ) - 45,000 [ krMJ'^-f^ ] 

, or 

(2-3) 

M, / (E,'B 2 *d 3 ) - 5 t m3.Mj-l-T-2 ] 


(2-4) 


For large magnets (E,> 3 uj). the mass of the force contaisent structnre, 
exceeds that of the condnotor winding, and 

- *»500 [ ] , or (2-5) 

(M^ + M 3 ) / (E.-B^-dg ) - 1 [ m3-MJ~l*f2 ] . ( 2 - 6 ) 

For a mostlj alHinna strnetnre, 63 * 3000. and a field of 5 T these equa¬ 
tions predict specific Msses of 600 and 120 kg/MJ respectiwely. 

It IS worth noting that the miniaom mass of the strnetnre according to 
Lewy's^ interpretation of the wirial theoreB cones to 11.5 kg/MJ for Amco 
21-6-9 stressed to 200,000 psi. Label 1 also graphs the specific masses for 
a nimiber of coils between 0.05 and 20 MJ and shows walnes between 600 and 
200 kg/MJ with a clear downtrend in the transition region. Fig. 1. 

For space rated OJ to 1 MJ 5 T aloiinna alloy solenoidal coils at 13 
K, I wonld thns suggest a conserwatiwe design valne of 

(M, + M 3 ) / E, - 300 kg/MJ . (2-7) 


D.2.3 Snrface Area of Inductor 


The Xol lowing is not an inductor design by any means, but merely some 
rough assoiptioms to arriwe at a plausible geoaetxy and a conservative upper 
1 iait tor the snrface area. They are: 

1. The inductor is a straight solenoid of rectangular cross-sect.on, 
with an inner radius a^. 

2. The length-to-diaaeter ratio ^ is 3. This approximates the 
energy density inside an infinite solenoid within 15%. The 
length is thus ^*s^, and the enclosed field voliae is 
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V = 6*jra^3 , (2-8) 

All the field energy is contained within that yolome so that the 
stored energy is 

E, * VB^ / 2-Hq . and (2-9) 

ai^ = Mo'Es f . (2-10) 

The average current density. J, in the winding volume is 100 
A/nuD^. or 10^ A/n^. so that the effective cross-sectional area of 
the conductor is 

Ao=I / J . (2-11) 

The umber of ampere turns per meter, n*I> are obtained from the 
long solenoid formula as 

n-I = B / pg . (2-12) 

The winding thickness is 

<il = n*I/J = B / (pg-J) . so that (2-13) 

a2 = *i'*‘‘*l ^l'*’ ® / (ho*^^ (2-14) 

a = 1 + B / (Pg-J'a^) . (2-15) 

There is a SO-naa space around the coil, so that 

Dewar Radius = 50 mm. (2-16) 

End supports and plumbing take up 300 mm, so that 

Dewar Length = iS-aj^+ 300 mm. (2-17) 

D-5 






D.2.4 LOSSES IN IHE INDDOOR 


ChMugiag fislds caass aaargj dissip«tioa is a. sap«rcoBdact or. Tho 
magnitude of these losses is strongly dependent on the detailed internal 
layout of the conductor wire and cable. Fig. 2 shows the 'nixed matrix* 
structure for a composite superconductor specifically designed for milli¬ 
second pulse applications. The superconducting filaments are made as fine 
as possible to reduce hysteresis lasses in the superconductor proper. Copper 
stabilizer is in close contact with each superconducting filament, but 
neighboring stabilizer elements are decoupled frcm each other by 'shields' 
of high resistiwity Cu-Ni alloy. Qis reduces eddy currents in the copper 
stabilizer and, in conjunction with a wire twist, limits similar induced 
'coupling' currents that are partly in the filaments and partly in the 
ma trix. 

The approximate relations from the prewious section can also be used to 
estimate the total conductor rolurne as 

- «*ra^*(4 • (2-a^ + d^) , (2-18) 

and that of the contained superconductor as 

Vg » . (2-19) 

Pinning centers are beliewed to enable superconductors to carry bulk 
supercurrents up to a critical current density For a wide range of 

fields the pinning force is constant, so that: 

Pinning Force =* x B . (2-20) 

Ihe hysteresis losses per unit wolume of superconductor represent the 
work done in pushing flux across pinning centers (Bean^), -Mid are for each 
half-cycle (charge or discharge); 

Q^/Vg - dB-Jg'dg / 4 [ J / ] , (2-21) 

where d i< the thickness of superconductor perpendicular to the changing 
field or the diameter of the superconducting filaments. 
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inproyed current densities ss well. 

Hie rate dependent coupling losses in mixed matrix condnctor that has 
been given the maximum practical twist pitch of eight wire diameters 
(Laquer*) are, 

®R/(Vo) = ( 3n/8)-*^*AB-Bdot/(p3-D) , (2-22) 

where: 

a = Wire diameter, typically 0.5 mm 
Bdot = dB/dt , Rate of field change [ T/$ ] 

P3 = p of high resistivity web, typically 10“^Q'm 
D = geometry dependent coupling parameter (0.148 max.). 

In applying these equations we shall assume that the entire volume of 
superconductor is subject to the full AB and then divide by two in Eq (2-21) 
and by four in Eq (2-22). A more detailed analysis, again requiring a more 
complete coil design, would include integration of the losses over the 
actual field distribution. 

D.2.5 Limitations on Fast Pulses 

The above loss equations are probably valid for field change rates of 
up to several thousand Tesla per second. For faster rates ai>d shorter 
pulses, however, one has to consider the distributed capacitance of the coil 
winding and the fact that the flux pinning mechanism, or its description 
breaks down for very high Bdot values. If the repetition rates are not too 
high it might be possible to extract longer pulses and then shape them in a 
non-superconducting device. 
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D.3.1 G«neral Design Priaciolds 


Hi* fuetioa of oryogoaic leads is to coadact s carreat froa s pover 
sooros St asbieat toapsrstare to sa eleetricsl dsrics st cxyogeaic teapers- 
tars. la doiag so tlisy also eoodaet boat beoroea tbess teapsratares. aad 
geasrata lossss. Siaos for a astal, tbexasl aad electrical coadactiTi- 

ties are faadaaeatal I7 related, ss sbova by tbs asar coastaacy of tbe Lorenz 
ratio, L-p k/T. there is ao way to bare oas witboat tbe other. 

The accepted procedare is to fabricate leads as ' eoaaterf 1 ow' beat 
ezcbaagezs aad size tbea so as to aiaiaize tbs beat flov iato tbs cryogenic 
eariroaieat at tbe rated cnrrsaC If tbe cryoelsctrical derics operates ia 
boiling liquid, all easrgy influx is taksa up ia tbe latsat beat of vapori- 
zatioa of the cryogea. All of the boil-off can bo channeled through tbe 
leads, and. proridsd tbe leads are tbezaal ly anchored in tbe cryogen, lead 
boil-off aad flov can be sel f-re gala ting. No beat will be conducted 
directly to tbe device. 

Miniana energy input froa veil-designed counterflcv leads into boiling 
belins saounts to 1.2 vatt/kA/lead, as sbovn eapirically aad froa nnaerous 
analyses, e.g., Ceilia aad Cl iaeako^ aad Efferson.^*^ 

D.3.2 Lead Refrigeration 

In an open systea vitb boiling liquid, no use can be aade of the 
sensible beat of tbe cold escaping gas. Most of tbe lead refrigeration is 
'free*. Hovever. ia a closed systea tbe sensible beat taken up by tbe lead 
flov gas ia returning to sabient teaperature also has to be prodnced by tbe 
refrigerator. The aore gas is 'bypassed* iato tbe carreat leads, froa tbe 
noaal return flov through tbe beat exchanger, tbe acre difficult it beccaes 
to aaintain refrigerator efficiency. 

Gerbold^^ has discussed tbe problea in general terns. Using bis anal¬ 
ysis together vitb seas unpublished Los Alaaosvork, I suggest a provisional 
foraula tor tbe coapressor pover, needed to refrigerate a pair of leads 

at a steady current 

Pli - 20*1 / Tg [ MCS ] (3-1) 
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wiier* ig IS tbs uiaimim gss toiiMrstBre in tne systea. It sbould be noted 
tbst this foxanls essnaos. in effect, 22% of Csrnot efficiency. It is 
probably optimistic end needs to be examined in more detail. Neyerthel ess, 
it serres to point out the treaendons refrigeration penalty imposed by high- 
current leads. A pair of steady state 50 kA leads at 12 C would require 
83 kW of ccapressor power. 

For quickly warying and pulse currents, it should be possible to time- 
awerage the current by relying on the heat capacity of the lead. However, 
depending on pulse leugth and repetition rate, such a scheme could introduce 
flow and pressure instabilities and degrade the heat exchange efficiency of 
the lead. 

D.3.3 Lead Resistance 

Current leads optimized for minimv refrigeration have a resistance 
that produces a fixed voltage drop of about 45 mV at the design current. 
This optimoB is quite independent of the material, its resistivity, the 
taiperature coefficient of resistivity, or the length of the lead. Nor does 
it matter whether the lead is of constant cross-section or tapered tcaard 
the cold end. We thus have for the optimaa resistance of a pair of leads: 

®L .opt “ 0-09/1 I MXS ] . (3-2) 

and for the corresponding power dissipation: 

Pl2 - 0.09-1 [ MIS ] . (3-3) 

A 50 kA current thus will dissipate 4.5 kW in a pair of leads. 

D.3.4 Lead Mass 

The length of each lead is usually set in the range of 0.5 to 1 m by 
space, mechanical, and insulation considerations. Since R«p*X /A. we get 
from eq'n (3-2) for each lead 

(>-/A) 3 0.045 / (I-P... ) [ MIS ] . (3-4) 







Claaxl7> » a«t«ri»l with th« lowast awarmgo rasistiwity. owar tha 

ralawant taapaxatnxa range, t.a.. a high parity natal such as A1 or Ca. 
woald than pamit tha lightest leads. Unfortunately, such leads are lass 
stable and lass tolerant of azcass cnrrents (nore likely to born oat) than 
those aada tron less pore natarials with higher resistiwities and saallar 
teaparstara coefficients. Iha increased mass and heat capacity of alloy 
leads also inprowes their theraal inertia for poise serrice, albeit at the 
expense of increased power dissipation when the cnrrsnt is increased above 
the optinoa. 

An aztreaa choice woald be brass with a resistivity of 8‘l0~^:*n at 
300 C and half that valaa at 10 L If the lead is tapered for a constant 
p/A ratio along its length. X. we get an average cross-sectional area 

Aj^g • 10-p^^g-X-I - C MB ] . (3-5) 

For 0.5 n. 50 kA brass leads the average crosssection woald be 300 ca^ and 
the mass of two leads woald be a nassive 261 kg. 

A better aatsrial woald be Deoxidized High Phosphoros Copper (DHP) with 
a resistance ratio of 7 and an average resistivity of l*10~^:*a. A pair of 
DHP 50 kA leads would weigh 27 kg if straight with a 30 cm^ cross-section, 
and 32 kg with a linear 5:1 taper froa 12 to 60 ca^. 

Inforaation for alaainaa alloy leads has not yet been collected, bat 
siailar resistivities are available. I therefore suggest an interia fosaula 
of 

Ml - 0.2*1 [ kg / kA ] . (3-6) 


Dewar vessels are passive stractores whose only purpose is to reduce 
the heat flow froa aabiont teaperatore to a cryogenic device, ezperiaent, or 
fluid. In the vacnna of space two nodes of heat transfer that are iaportant 
in an ataospheric environaent are absent. These are gas convection and gas 
conduction. Only radiation heat transfer between surfaces and condnction 
through solids reaain. 

The rate at which a surface of area, A. eaits therasl radiation is 
given by the Stefan^Bol tzaann relation 
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P - a*«’A*T* 


(4-1) 


whero: s is the taissiTity of the surface, and 

a - 5.67*10“* [ W*«“^*r“^ ] the radiation constant. 

The net exchange of radiant energy between two snrfaces of equal (low) 
enissiTity and equal area at and talces place at 

Pj- o*s*A*(Ti^ - Tq^) . (4-2) 

Clearly, radiative heat transfer is controlled alnost entirely by the 
high teaperatnre surface. A 300 C alvinia surface with e-0.02 will radiate 
9.2 W/b 2. The r eturn radiation froM a surface at 10 or 20 E is lower by 
four or five orders of nagnitude, and therefore is insignificant in reducing 
the heat leak. 

Conduction through solids, on the other hand, is expressed by the 
Fourier equation, which is linear in the teaperature gradient: 

Pj- k*A*(T^ - T^) / X . (4-3) 

Optiaally designed dewars usually have equal heat losses through radia- 
tiou and conduction. There is no sense in trying to make one heat leak 
orders of magnitude lover than the other. As a matter of fact, any signifi¬ 
cant additional energy losses in the systea should also be considered when 
designing the dewar insulation. 

D.4.1 Multilayer Insulation 

The ideal insul ation vould consist of a niaber, N, of concentric, mass¬ 
less, loveaissivity, isothexaal. floating radiation shields in a perfect 
vacuua. Such an arrangeaent would reduce the heat leak given by equation 
(4-2) by a factor of N+1. 

Coaaercial multilayer insulation, known as 'super-insul ation', tries to 
approximate the ideal by interleaving thin aliainna foils with fiberglass 
paper spacers or by using multiple wraps of crinkled alosinized mylar. The 
latter is the lighter option and also can be aliainixed with strips or small 
squares so as to minimize eddy currents induced by changing magnetic fields. 
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Hia teaperature profile ia these iasolstioas. howerex, is aot ss oae sight 
e^eot tor s series of pare rsdistioa shields, bat becoaes 1 iaear below 70 
X. ^is iadicates that coadactioa throagh eoataots is predoainaat ia the 
low teaperatare regioa. effectiToly shortiag oat the iasalatioa that woald 
be proTided bj the radlatioa shields. Oa the other haad, with the abseace of 
grawity ia a space dewar, it is eoaceirable that coatact forces aad the 
resaltiag heat flow coald be reduced throagh proper desiga. 

Altaiaized a^lar iasalatioa is awailabls ia balh fora. Applicatioa 
data, pablished by Natioaal Research Corpora tioa^^, saggest that the 
aaterial is aost effectiwe whea applied at 16 layers/ca (40 layers/iach) aad 
thea will hare a heat flax be tweea 300 E aad 20 E of 

?! - 18.9 / N [ ra“^*layors"^ ] . (4-4) 

prowided the aiaber of l^ers, N, is be tweea 30 aad 150. 

The aass of the iasalatioa is 

Mj - 8.94'10”^*N [kr»"^*l»y«rs ] . (4-5) 

As loag as the dewar will oaly be operated ia space, its iaaer al laiiana 
tl 1 oy she! 1 aeed be ao aore thaa 5 oa thick aad its oater shel 1 a mere 1 nna 
This woald gire a aass of 


Md - 16J I kg’a-2 ] . (4-6) 

For short palses, eddy carreats ia the dewar shell aay be aadesirable aad a 
fiberglass epoxy dewar of siailar aass aay haw e to be used. 


laforaatioa oa space rated refrigerators is rather scarce ia the opea 
litaratare, probably for classifiertioa aad proprietary reasoas. However, 
Strobridge^^ has dooa a thorough survey of cryogeaic refrigerators. After 
takiag out the inevitable (]arnot factor 

”0 - / (T. - T^) . (5-1) 
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the otabiaed theoael end aechanicel efficiency, does not depend on the 

teaperetnze, T^, between 1,8 end 90 K. Fig. 3 telwn fran his paper shows 
that it is only a fnnction of the refrigeration capacity (frigiwatts), Pq. 
For saall refrigerators with large snrface-to-wolone ratios is saal 1. It 
reaches IS % for kilowatt aachines and levels off at about 33 % for aegawatt 
units. 

Strobridge also gives a relation for costs, (which have to be 

adjusted for inflation since 1969) as a fnnction of the compressor power at 
aabient teaperature, 


Cj - dS'P^O*"^ (1969 i, watt ] . (5-3) 

Finally, we can extract the following relations frca his plots of aass and 
volnae against frigiwatts for refrigerators in the 10 to 30 E range: 


Mr - 

gO-p^O.T ( kg. 

watt ] , and 

(5-4) 

Vr * 

0.0SS-Pg®-*^ t 

B^, watt ] . 

(5-5) 


D.5.1 Refrigeration Alternatives 

Stored cryogens present an alternative to aechanicel refrigeration, but 
only for relatively short aissions and low refrigeration powers. The latent 
heat of vaporization can be used effectively above the triple point tempera¬ 
ture, Tp, Various cryogens have been used in this wanner for cooling either 
equipment or radiation shields. For liquid hydrogen the cooling power 
amounts to about 4S0 kJ/kg. The triple point is at 13.8 C, mandating a 
slight increase in the operating temperature of the inductor or a short span 
(aagnetic 7) refrigerator. Detailed analysis of these options is beyond the 
scope of this note. 


We soamarize the use of the parmaetric model with the example of a SOO 


kj, SO kA. 5 T superconducting 


:orage indui 


at 13 E and refrigerati 








12 I. rocogaiziag that aost of tha aaabors are iateaded to be coaserrative. 
i.e. . oa the high side. 

Froa 2-7 we get the aass 

Mt + Ms - 150 kg. 

Froa 2-10 laaer radiaa: aj ■ 0.139 a 

Coil leagth: 2‘b - ^‘a^ - 0.832 a 

Froa 2-13 Wiad. thickas: d]^ * 0.040 a 

Froa 2-14 Oater radiaa: * 0.139 0.040 >.179 a 

Froa 2-16 Dewar radios: a2'^ 0-05 ■ 0.229 a 

Froa 2-17 Dewar leagth: 2*b 0.3 > 1.132 a 

Dewar sarface: Aq “ 1*W 

Froa 2-18 Coad. Voltaie: > 3.33*10“^ a^ 

This does not iaclode aay straotare aad correspoads to a aass of sboot 
100 kg coapared to the abowe coaserwatiwe ▼aloe of ISO kg. 

Froa 2-19 Sapercoad Voloae 

Vq - 3.33-10"3 a3 

Froa 2-21 Hysteresis Loss for AB > S T 
do • 5 iui 

Qq - 21 J / half cycle 

if the eatire coadactor sees the fall field chaage. otherwise divide by 2 
aad get 10.S J. 

At a rep rate of 1 Hz this dissipates 

Pg » 2’IO.S > 21 frigiwatts. 

(Saae loss for charge aad discharge). 

Froa 2-22 Coapl iag Loss for the givea typical valoes aad 
AB - 5 T 

Bdot > SOOO X/s (1 as palse) 

> 1660 J / discharge cycle 

if the eatire coadactor sees the fall field chaage, otherwise divide by 4 
aad get 415 J. 

At a rep rate of 1 Hz this dissipates 

Pj > 415 frigiwatts. 

(Assniag that Bdot for recharge is negligible). 
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From 3-1 Lead Refrigeration for SO kA 

Pli “ 83 kW Compressor Power 
Divide by 2 for average current of 25 kA. 

PLi - 42 kW 

From 3-3 Lead Dissipation for 50 kA 

Pl 2 “ M of Pnlse Power 
Divide by 4 for average power dissipation. 

From 3-6 Lead Mass for SO kA 

Ml - 10 kg 

Soi of average low temperature power dissipation from hysteresis and 
coupling (neglecting 'ordinary* eddy currents and deferring tbe lead prob¬ 
lem) 

Pq+Pl* 436 frigiwatts 

If these losses are present all the time, dewar losses would need be no 
less than 44 w. 

From 4-4 Layers of Snperinsnlation to achieve 44 W 

N =• 18.9'A0/44 - 0.84 , 

or a single radiation shield. 

A simple 40 layer insulation system would have a heat leak of 
Pj= 18.9*Ao/40 =0.93 frigiwatts. 

This would minimize standby power requirements. 

From 4-5 and 4-6 Dewar Mass 

Ml+Mp - 0.7+31.8 = 32.5 kg. 

From 5-1 Carnot Factor at 12 K 

nc = 0.042 or 4.2'!b 

From Fig. 3 Refrigerator Efficiency (iterated) 

** 0.13 or 13% 

From 5-2 Compressor Power (neglecting leads) 

P,= 437/(0.042-0.13)= 80 kW 

Total Compressor Power 

Pli 122 KW 

From 5-2 and 5-4 combined Refrigerator Mass 
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Total Mass 


5680 kg 


Mr- 150+10+32+5680 - 5872 kg 

To calculato efficiaaciss sod i |2 Boxa infonaation is oaadad on oparation 
and dnty cyclas and on tka priaia povar supply. 

Claarly. the biggest power and nass reqaireaants are associated with 
the size of the refrigerator, which, in torn, is determined by the currents 
in the leads and the coupling (eddy current) losses in the composite supez- 
conductor. 
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0.8 LIST OF FI6DKES 

Fig. 1. Specific Mass of Sapercondncting Coils, from Lnbell (1973) 

Fig. 2. Mixed Matrix Super conduct or Designed for Millisecond Pulse Applica¬ 
tions. from Laquer (1976) 

Fig. 3. Efficiency of Low Temperature Refrigerators as a Function of 
Refrigeration Capacity (Frigiwatts). from Strobridge (1969) 
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